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ABSTRACT 

Observational  evidence  of  seasonal  variability  below  the  main  thermo- 
cline  in  the  eastern  North  Atlantic  is  described,  and  a  theoretical  model 
of  oceanic  response  to  seasonally  varying  windstress  forcing  is  constructed 
to  assist  in  the  interpretation  of  the  observations.  The  observations  are 
historical  conductivity-temperature-depth  data  from  the  Bay  of  Biscay 
region  (2°  to  20°W,  42°  to  52°N) ,  a  series  of  eleven  cruises  over  the  three 
years  1972  through  1974,  spaced  approximately  three  months  apart.  The 
analysis  of  the  observations  utilizes  a  new  technique  for  identifying  the 
adiabatically  leveled  density  field  corresponding  to  the  observed  density 
field.  The  distribution  of  salinity  anomaly  along  the  leveled  surfaces  is 
examined,  as  are  the  vertical  displacements  of  observed  density  surfaces 
from  the  leveled  reference  surfaces,  and  the  available  potential  energy. 
Seasonal  variations  in  salinity  anomaly  and  vertical  displacement  occur 
as  westward  propagating  disturbances  with  zonal  wavelength  390  (±50)  km, 
phase  71  (±30)  days  from  1  January,  and  maximum  amplitudes  of  ±30  ppm  and 
±20  db  respectively.  The  leveled  density  field  varies  seasonally  with  an 
amplitude  corresponding  to  a  thermocline  displacement  of  ±15  db. 

The  observations  are  consistent  with  the  predictions  of  a  model  in 
which  an  ocean  of  variable  stratification  with  a  surface  mixed  layer  and 
an  eastern  boundary  is  forced  by  seasonal  changes  in  a  sinusoidal  wind- 
stress  pattern,  when  windstress  parameters  calculated  from  the  observa¬ 
tions  of  Bunker  and  Worthington  (1976)  are  applied. 
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Chapter  1 

1.1  Introduction 

Seasonal  variations  in  the  strength  of  the  wind  are  comparable  in 
magnitude  to  the  mean  wind  over  much  of  the  world's  ocean.  Consequently, 
significant  seasonal  fluctuations  in  the  wind  driven  ocean  circulation 
might  be  anticipated.  The  ocean  adjusts  to  periodic  forcing  at  the  sur¬ 
face  through  a  combination  of  barotropic  and  baroclinic  motion.  The  baro- 
tropic  motion,  which  is  independent  of  the  stratification,  represents  the 
response  of  the  water  column  as  a  whole  to  the  changes  in  the  surface 
mass  field  which  occur  as  a  result  of  the  stress  exerted  on  the  surface 
by  the  wind.  Baroclinic  motion  exists  only  if  the  fluid  is  stratified, 
and  represents  the  adjustment  of  the  density  field  to  the  imposed  forcing. 
The  time  scale  for  baroclinic  adjustment  is  longer  than  for  barotropic 
adjustment.  Veronis  and  Stommel  (1956)  demonstrated  that  the  seasonal 
time  scale  represents  a  crossover  point  between  high  frequency  forcing 
with  dominantly  barotropic  response  ana  very  low  frequency  forcing  with 
dominantly  baroclinic  response.  For  seasona)  forcing,  the  barotropic 
and  baroclinic  components  are  of  approximately  equal  importance,  so  that 
the  adjustment  of  the  ocean  to  seasonally  varying  windstress  is  not 
confined  to  the  surface  layers.  In  addition  to  seasonal  variations  in 
the  wind,  there  is  a  seasonal  cycle  to  the  heating  and  cooling  of  the 
ocean  surface,  which  may  cause  significant  seasonal  variations  in  the 
large  scale  thermohaline  circulation  of  the  ocean. 

The  general  context  of  this  work  is  the  presentation  and  inter¬ 
pretation  of  observational  and  theoretical  evidence  for  seasonal  vari¬ 
ability  in  the  ocean  at  depths  below  the  region  of  the  direct  atmos.- . 
pheric  influence.  The  observational  evidence  i'svtaken  from  a  three  year 


8 


series  of  hydrographic  cruises  in  the  eastern  North  Atlantic.  A  theoret¬ 
ical  model  of  the  oceanic  response  to  seasonal  windstress  forcing  is 
constructed  using  the  observed  oceanic  parameters,  and  the  model  pre¬ 
diction  is  compared  with  the  observed  variability. 

Veronis  and  Stommel  (1956)  demonstrated  that  seasonal  forcing 
by  a  mov.ng  wind  stress  pattern  in  an  unbounded  ocean  results  in  very 
small  vertical  and  horizontal  displacements  in  the  main  thermocline.  In 
that  case  the  length  scale  of  the  forcing  and  response  are  the  same. 
However,  the  presence  of  lateral  boundaries  imposes  an  additional  con¬ 
straint  cn  the  system  of  no  flow  into  the  boundary.  In  general,  the 
forced  response  of  the  unbounded  ocean  cannot  satisfy  that  boundary  con¬ 
dition,  so  that  free  motions  (solutions  to  the  unforced  equations)  must 
be  added  to  the  forced  solution.  If  the  system  of  equations  is  linear, 
the  superposition  of  free  and  forced  motions  which  satisfies  the  lateral 
boundary  condition  is  also  a  solution  to  the  governing  equations.  The 
scales  of  the  free  motions  are  determined  by  the  geometry,  the  frequency, 
the  stratification,  and  the  lateral  boundary  condition.  In  certain  cases, 
the  dominant  free  response,  which  necessarily  has  a  wave  amplitude  com¬ 
parable  to  that  of  the  forced  response,  may  have  a  horizontal  scale  much 
smaller  than  that  of  the  forcing.  In  that  circumstance,  the  vertical  and 
horizontal  particle  displ' cements  of  the  dominant  free  response  will  be 
larger  than  those  of  the  forced  response.  This  is  due  to  the  geometrical 
fact  that  the  same  amount  of  en  rgy  (which  is  proportional  to  the  square 


of  the  wave  amplitude)  distribute.':  over  a  smaller  area  will  result  in 


larger  displacements.  If  the  free  response  scale  is  small  enough,  the 
response  of  the  ocean  to  large  scale  atmospheric  forcing  can  be  detected 
observationally,  even  though  the  forced  response  alone  is  too  weak  to  be 


observed. 


The  primary  hypothesis  of  this  work  is  that  the  observed  annual 
variations  in  wind  stress  curl  over  the  North  Atlantic,  in  combination 
with  a  meridional  barrier  at  the  eastern  edge  of  the  ocean,  are  capable 
of  producing  seasonal  variability  below  the  main  thermocline  of  suf¬ 
ficient  amplitude  to  be  detected  observationally  in  the  eastern  part  of 
the  basin.  There  is  a  corollary  hypothesis  that  low  frequency  signals 
wnich  involve  small  particle  velocities  can  be  observed  indirectly  if 
the  motion  occurs  in  a  region  of  strong  horizontal  gradient  of  any  tracer, 
by  following  the  movement  of  parcels  tagged  by  an  initial  tracer  con¬ 
dition.  The  strong  ambient  gradient  acts  to  simplify  the  signal,  in  the 
sense  that  lateral  displacements  will  result  in  anomalies  from  the 
initial  state  of  the  tracer  involved.  Therefore,  information  from  tracer 
distributions  along  density  surfaces  can  be  used  to  detect  velocity  sig¬ 
nals  too  small  to  be  resolved  by  standard  Eulerian  current  measurements, 
provided  the  lateral  tracer  gradients  are  sufficiently  strong. 

Conceptually,  this  work  is  concerned  with  the  transfer  of  energy 
from  large  scales  in  the  atmosphere  to  considerably  smaller  scales  in  the 
ocean.  Although  the  details  of  the  transfer  process  itself  are  not 
studied,  the  estimation  of  the  observed  available  potential  energy  in  the 
ocean  is  discussed  at  length.  The  available  potential  energy  or  APE  is 
that  part  of  the  total  potential  energy  which  is  actually  available  for 
conversion  to  kinetic  energy.  Formally,  it  is  defined  as  the  difference 
between  the  horizontally  averaged  total  potential  energy  and  the  adiabatic 
minimum  in  potential  energy  obtained  when  surfaces  of  constant  potential 
density  coincide  with  geopotential  surfaces  (i.e.  potential  energy  cor¬ 
responding  to  the  adiabatically  leveled  reference  density  field) .  Pre¬ 
vious  estimates  of  APE  in  which  effects  of  compressibility  are  ignored 
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and  reference  states  other  than  the  adiabatically  leveled  field  are  used 
do  not  allow  the  precise  evaluation  of  errors  in  the  estimate. 

A  computational  technique  for  determining  the  adiabatically  leveled 
reference  density  field  using  the  full  equation  of  state  is  described 
below  whose  accuracy  is  limited  only  by  the  sampling  and  measurement 
errors  of  the  observations.  The  technique  permits  the  evaluation  of  the 
effects  of  measurement  and  finestructure  errors  on  the  calculation  of  APE, 
as  well  as  the  errors  which  result  from  neglecting  compressibility  in  the 
Boussinesq  estimate  of  APE.  Because  geopotential  is  not  measured  accur¬ 
ately  in  the  ocean,  pressure  surfaces  rather  than  geopotential  surfaces 
are  used  in  the  definition  of  the  reference  state.  The  leveled  density 
surfaces  are  also  the  correct  surfaces  for  examining  tracer  distributions, 
and  are  employed  in  all  calculations  involving  tracers  in  this  work.  The 
adiabatic  leveling  technique  provides  a  precise  and  consistent  overall 
framework  for  the  analysis  of  hydrographic  data  in  terms  of  both  energetics 
and  tracer  distributions. 


1.2  Earlier  Work 


There  is  little  observational  evidence  for  seasonal  variability 
in  the  oceans  below  the  surface  layers.  The  reason  for  this  is  two-fold: 
only  recently  have  observations  of  sufficiently  long  duration  been 
available,  and,  except  at  low  latitudes,  the  expected  annual  signal  is 
small  compared  to  the  erergetic  western  boundary  currents  and  mesoscale 
eddies.  In  the  near  equatorial  region  where  the  baroclinic  response 
time  of  the  ocean  is  shorter  than  at  higher  latitudes,  seasonal  variability 
has  been  observed.  White  (1977)  found  evidence  for  propagation  of  long 
baroclinic  Rossby  waves  in  the  main  thermocline  (depth  about  200  m)  of 
annual  period  in  MBT  (mechanical  bathythermograph)  data  from  the  tropical 
North  Pacific.  From  the  phase  information  he  inferred  that  the  source  of 
the  waves  was  the  eastern  boundary.  In  a  similar  analysis  White  (1978) 
presented  evidence  from  the  mid-latitude  North  Pacific  for  seasonal  fluc- 
tuations  in  the  depth  of  the  main  thermocline.  He  demonstrated  that  the 
phase  of  those  fluctuations  matched  the  phase  of  the  observed  windstress 
curl;  however  the  observed  amplitude  was  5  to  10  times  that  expected  from 
the  theory  of  Veronis  and  Stommel  (1956) .  No  explanation  of  the  amplitude 
mismatch  is  given  by  White. 

In  the  North  Atlantic  the  main  thermocline  is  considerably  deeper 
than  200  m  in  mid-latitudes  (600-800  m) ,  well  removed  from  any  direct 
solar  influence.  There  are  also  major  sources  of  deep  and  bottom  water 
whose  formation  at  the  surface  at  high  latitudes  and  in  the  Mediterranean 
Sea  occurs  seasonally.  In  the  western  part  of  the  basin  there  is  an 
energetic  mesoscale  eddy  field  associated  with  the  Gulf  Stream  system. 
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In  the  topographically  confined  Florida  Current  whose  transport  con¬ 
tributes  significantly  to  the  Gulf  Stream,  Niiler  and  Richardson  (1973) 
cite  evidence  of  seasonal  fluctuations  representing  half  the  total  varia¬ 
bility  of  the  current  and  about  10%  of  its  total  transport.  Away  from  the 
Florida  Strait,  however,  no  evidence  for  seasonal  variability  below  the 
main  thermocline  has  been  cited  in  the  western  North  Atlantic.  The 
energetic  eddy  field  associated  with  the  Gulf  Stream  will  tend  to  mask 
the  relatively  weak  annual  signal  expected  in  the  interior,  necessitating 
long  observational  records  to  detect  the  signal  unambiguously.  Wunsch 
(1972)  computed  spectra  of  temperature  and  dynamic  height  anomaly  using 
the  13  years  of  PANULIRUS  hydrographic  data  taken  semi-monthly  at  a  single 
location  near  Bermuda.  He  found  a  pronounced  annual  peak  in  temperature  at 
10  m  depth,  a  less  pronounced  peak  at  100  m  and  none  at  800  m.  He  also 
found  annual  peaks  in  the  dynamic  height  anomaly  spectra,  but  all  the 
energy  at  that  frequency  comes  from  depths  shallower  than  200  m.  Thus, 
even  in  a  very  long  single  record  in  the  western  North  Atlantic  the  sea¬ 
sonal  signal,  if  it  exists,  has  not  been  extracted  from  the  noise. 

Nevertheless,  the  available  meteorological  observations  in  the 
North  Atlantic  indicate  that  there  is  a  strong  annual  cycle  to  the  wind- 
stress  curl,  with  large  horizontal  scales  (Bunker  and  Worthington,  1976) 
so  that  some  annual  response  in  the  ocean  interior  is  expected.  It  is 
anticipated  that  the  observation  of  that  annual  signal  is  most  likely 
to  be  made  in  the  eastern  basin,  away  from  the  energetic  Gulf  Stream 
system,  and  particularly  in  a  region  of  relatively  weak  stratification 
near  the  eastern  boundary.  As  discussed  briefly  in  the  introduct ' ^ , 
the  boundary  results  in  free  waves  as  part  of  the  solution  to  the  forced 
problem;  the  weaker  the  stratification,  the  shorter  the  horizontal  scale 


of  the  free  waves,  which  in  turn  results  in  larger  displacements. 

Bryan  and  Ripa  (1978)  constructed  a  model  of  the  oceanic  response 
to  large  scale  continuous  wind  forcing  at  low  frequencies,  for  the 
special  case  of  a  flat  bottomed  ocean  with  depth  variable  stratification 
and  a  single  meridional  barrier  at  the  eastern  edge.  They  applied  their 
model  to  the  mid-latitude  North  Pacific  at  frequencies  corresponding  to  periods  of 
3  and  6  years,  using  an  idealization  of  the  observed  windstress,  and 
appropriate  oceanic  parameters.  They  calculated  the  resultant  scales 
of  the  dominant  free  modes  and  estimated  the  apparent  vertical  pro¬ 
pagation  for  comparison  with  observational  estimates  of  the  vertical 
propagation  of  low  frequency  temperature  anomaly  structures  in  the 
North  Pacific.  Bryan  and  Ripa  made  no  attempt  to  calculate  the  amplitude 
of  the  response,  although  they  discuss  the  phase  of  the  solution  at 
length.  Their  approach  is  appropriate  to  the  present  work,  and  their 

model  is  discussed  in  detail  when  it  is  applied,  with  some  modification 

0 

and  extension,  to  the  three  years  of  CTD  data  from  the  eastern  North 
Atlantic,  in  Chapter  3. 

The  region  of  the  eastern  North  Atlantic  in  which  the  data  used 
in  this  work  were  collected  has  no  major  sub-surface  current  systems, 
with  the  possible  exception  of  a  weak  (.01  to  .05  m  sec  *)  poleward 
eastern  boundary  current  of  60  to  250  km  width  and  undetermined  depth 
range.  Away  from  the  coast  the  available  direct  current  measurements  in¬ 
dicate  very  low  mean  flows  (^  .01  m  sec  1)  (Swallow  et  al.,  1977).  There 
is  no  evidence  for  a  strongly  energetic  eddy  field  as  is  found  in  the 
western  basin,  although  isolated  eddy-like  features  have  been  observed 
(Swallow,  1969  and  Gascard,  1980)  with  velocities  of  .10-. 20  m  sec  1 
associated  with  the  cyclonic  flow,  and  some  indication  of  slow  westward 


14 


drift,  on  the  order  of  .02  m  sec  1 . 

There  are  a  number  of  water  masses,  as  distinguished  by  temperature- 
salinity  relationship,  found  in  this  region.  Of  primary  interest  to  this 
work  are  the  water  masses  found  between  the  main  thermocline  and  roughly 
2000  m,  which  is  the  deepest  that  the  CTD  stations  penetrate.  Those  water 
masses  are  Mediterranean  Outflow  Water  (MOV?) ,  which  is  found  in  the 
western  Mediterranean  Sea  and  is  characterized  by  high  salinity  as  a 
function  of  temperature  compared  to  North  Atlantic  Deep  Water,  and 
Labrador  Sea  Water,  formed  in  the  Labrador  Basin  and  similarly  charac¬ 
terized  by  low  salinity.  The  interaction  of  these  two  water  masses 
within  the  area  studied  results  in  strong  mean  lateral  gradients  of 
salinity  anomaly  (defined  as  the  difference  between  the  observed  salinity 
and  a  reference  salinity) .  Other  work  has  been  done  which  exploits 
these  strong  lateral  (i.e.,  isopycnal)  gradients  of  salinity  anomaly, 
notably  Needier  and  Heath  (1975)  and  Katz  (1970).  Of  principal  interest 
here  is  the  identification  in  both  cases  of  strong  lateral  gradients  of 
a  suitably  defined  salinity  anomaly. 

In  order  to  study  the  lateral  distribution  of  any  tracer  (such 
as  salinity  anomaly)  it  is  necessary  to  define  the  density  surfaces  to 
be  used.  Density  surfaces  rather  than  horizontal  surfaces  are  used 
since  a  minimum  of  work  is  performed  in  moving  a  parcel  along  a  surface 
of  constant  density.  Ideally  what  is  sought  are  those  level  surfaces 
which  are  connected  to  the  observed  density  field  by  adiabatic  displace¬ 
ments  at  each  depth.  Only  recently  has  a  computational  technique  been 
developed  for  determining  these  implicit  adiabatically  leveled  surfaces 
(Bray  and  Fofonoff ,  1980) .  That  technique  is  described  in  section  .4 


of  this  chapter. 
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Montgomery  (1938)  in  his  isentropic  analysis  of  the  southern 
North  Atlantic  used  0  surfaces  as  an  approximation  to  constant 
potential  density  surfaces.  (If  a  parcel  of  water  it  pressure  p  has 
temperature  T  and  salinity  S,  then  0  is  giver  by 


0  (p,T,S)  =  (p(p  ,T,S)  -  1)  X  103 

tl  cl 


p  the  density,  evaluated  at  atmospheric  pressure  p  .)  Surfaces  of 

cl 

constant  a  are  a  reasonably  good  approximaton  to  constant  potential 
density  surfaces  for  observations  shallower  than  10C9  m.  A  better  ap¬ 
proximation  is  a  surface  of  constant  0g,  for  which  the  effect  of  the 
adiabatic  change  of  temperature  with  pressure  has  been  corrected.  If 
6  (P/Pr)  is  the  potential  temperature  at  pressure  p  referred  to 
pressure  p^  (Fofonoff,  1977): 


9  (p»Pr)  =  t(p)  + 


r  dp 


(T  the  adiabatic  lapse  rate 
then  is  defined  by 


<!>*> 


0Q  (p,T,S)  =  (p(p  ,  0  (p,p  ),  S)  -  1)  x  103 
o  a  a 


p  the  atmospheric  pressure.  However,  in  the  deep  water,  0q  surfaces 
are  no  longer  a  good  approximation  to  constant  potential  density  surfaces 
(due  to  large  changes  in  the  coefficient  of  thermal  expansion  with 
pressure) ,  so  that  some  pressure  other  than  atmospheric  pressure  must 
be  used  as  the  reference.  Reid  and  Lynn  (1971)  used  0^,  0^  and  0^, 
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The  adiabatic  leveling  technique  uses  a  different  reference  pressure  p 
at  each  depth  for  which  a  surface  is  calculated,  with  the  additional 
constraint  that  mass  is  conserved  in  the  leveling  process.  The  questions 
of  the  appropriate  p  for  a  given  depth  range  and  of  transition  from 
one  Pr  to  another  are  thus  separated. 

The  adiabatically  leveled  reference  state  also  represents  the 
adiabatic  minimum  in  potential  energy:  that  is,  no  further  energy  can 
be  extracted  from  the  system  by  adiabatic  processes.  The  potential  energy 
in  the  observed  field  which  is  actually  available  for  conversion  to 
kinetic  energy  (available  potential  energy  or  APE)  is  the  difference  in 
the  horizontally  averaged  total  potential  energy  and  the  adiabatic  mini¬ 
mum  of  potential  energy.  Lorenz  (1955)  developed  this  concept  in  ap¬ 
plication  to  the  atmosphere.  If  the  atmosphere  is  taken  to  behave  as  an  ideal 
gas,  the  total  potential  energy  and  internal  energy  are  proportional; 
consequently  Lorenz  defines  APE  as  the  difference  in  enthalpy  (the  sum 
of  potential  and  internal  energies)  between  the  observed  and  reference 
states.  In  a  compressible  fluid  the  vertical  displacement  of  an  isopycnal 
may  be  accompanied  by  changes  in  volume  due  to  differential  compressibility. 
(Compressibility  is  a  function  of  pressure  and  temperature.)  In  the 
atmosphere  those  volume  changes  are  large  (one  third  of  the  total  enthalpy 
change)  and  act  to  reduce  the  potential  energy.  In  the  ocean  they  are 
small  and  generally  act  to  increase  the  potential  energy.  As  will  be 


demonstrated  in  some  detail  in  the  next  section,  these  changes  in 
potential  energy  result  in  very  small  changes  in  APE,  since  most  of  the 
change  is  not  available  for  conversion  to  kinetic  energy  at  all,  but 
must  go  into  changing  the  reference  field. 

A  difficulty  of  definition  then  arises  as  to  what  represents  the 
true  APE.  It  is  neither  the  enthalpy  difference  nor  the  change  in 
potential  energy  which  includes  all  of  the  change  due  to  compressibility 
In  the  following  section  a  system  of  nomenclature  is  proposed  to  distin¬ 
guish  the  different  estimates  of  APE;  a  simple  example  is  discussed 
which  illustrates  the  problem  in  more  detail,  and  the  computational  tech 
nique  for  identifying  the  adiabatically  leveled  reference  density  field 
is  presented. 


1 . 3  Potential  Energy 


(The  discussion  in  this  section  follows  the  similar  discussion  in 
Bray  and  Fofonoff ,  1980) . 

The  total  potential  energy  contained  in  a  column  of  seawater  of 
unit  cross-sectional  area  relative  to  a  reference  level  z is  given  by 

z 


E  = 


P(<f>  -  <j>r)  dz 


1.1 


z 


r 


with  p  the  density,  4>  the  geopotential  and  zg  the  height  of  the 
surface  of  the  ocean.  If  the  fluid  is  in  hydrostatic  balance,  (1.1) 
may  be  rewritten  in  terms  of  pressures 


since 


(<l>  -  <l>r>  dp 


t 

dp  =  -pg  dz  =  -p  dcj> 


The  geopotential  <J>  may  also  be  written  in  terms  of  pressures 


with  a  the  specific  volume,  a  =  1/p.  Thus,  in  terms  of  pressure  and 
specific  volume  the  total  potential  energy  is 


E  = 


a  dp "dp. 
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Lorenz  (1955) ,  in  his  studies  of  the  general  circulation  of  the 
atmosphere,  pointed  out  that  the  total  potential  energy  has  little 
dynamical  relevance,  since  only  a  very  small  percentage  is  available 
for  conversion  to  kinetic  energy;  specifically,  that  which  is  in 
excess  of  the  adiabatic  minimum.  That  excess,  for  a  column  of  unit 
cross-section  is 


Ae  = 


(cu  -  af)  dp"  dp 


ps  p 

with  af  the  reference  (adiabatically  leveled)  specific  volume,  a 
function  only  of  pressure,  and  the  observed  specific  volume.  AE 

can  be  either  negative  or  positive,  depending  upon  the  sign  of  ou  -  ctf- 
However  there  is  a  net  positive  storage  of  available  energy  averaged 
over  a  volume  for  which  mass  is  conserved  during  leveling.  The  net 
storage  is  the  total  available  gravitational  potential  energy  AE, 


averaged  horizontally  over  an  area  A: 


Ae  dA  = 


(ai  ”  af*  dp^  dp 


(1.2) 


A  Ps  P 

Contributions  to  the  total  AE  from  within  the  column  are  identified 
as  the  total  available  gravitational  potential  energy  (TGPE) : 


TGPE (p) 


(1.3) 


with  units  of  potential  energy  per  unit  area  (J  m  2)  and  contributions 
to  TGPE  at  each  pressure: 
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GPE(p)  =  g 


d  TGPE 


- 

-  =  (ct.  —  a  )  dp-* 
J  it 


(1-4) 


with  units  of  potential  energy  per  unit  mass  (J  kg  1 ) . 

As  an  illustration  of  the  relationship  between  the  gravitational 
potential  energy  and  the  thermodynamic  energies  of  the  system  consider 
the  following  simple  example.  A  rigorous  mathematical  derivation  is 
found  later  in  section  3.  Take  a  volume  of  seawater  of  length  and 
width  L  and  infinite  vertical  extent,  with  specific  volume  initially 
level.  Next,  displace  two  columns  of  fluid  of  unit  mass  adiabatically 
a  distance  {in  pressure  units)  tt  :  one  upward  and  one  downward,  there¬ 
by  conserving  mass  at  each  pressure.  If  the  reference  density  varies 
lineally  with  depth,  no  net  work  is  done  (except  small  changes  of 
internal  energy,  which  are  discussed  later) .  That  is,  the  horizontal 
average  of  "  af  is  zero  since  the  perturbation  of  the  upward 
displacement  results  in  the  same  change  of  volume  as  that  of  the 
downward  displacement,  unless  there  are  gradients  of  compressibility 
present,  in  which  case  conversions  of  internal  energy  to  potential 
energy  are  possible.  However,  if  density  is  a  non-linear  function  of 
depth,  the  upward  displacement  results  in  a  different  change  of  volume 
than  the  downward  displacement;  there  is  a  net  change  in  specific 
volume  averaged  over  L  ,  and  net  work  is  performed.  A  series  of 
schematics  relating  to  this  example  are  found  in  Pig.  1.1.  Although 
the  motion  is  presumed  to  extend  throughout  the  fluid  the  representation 
of  the  observed  and  reference  surfaces  as  functions  of  pressure  and 
specific  volume  are  shown  for  a  single  level.  In  Fig.  1.1a  the  level 


: 
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Figure  1.1  Schematic  of  displacements  in  an  initially  level  steric 


field.  See  text  for  detailed  explanation. 
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field  is  shown.  In  Fig.  1.1b  the  level  field  is  perturbed  by  displace¬ 
ments  ir.  The  amount  of  work  done  can  be  estimated  as  follows:  the 
change  in  volume  at  pressure  p^  for  the  upward  displaced  (— TT ) 
column  is  given  approximately  by 


Aa 

U 


(pf  +  tt) 


whereas  the  change  in  volume  corresponding  to  the  downward  displacement  is 


da 

dp 


-  it) 


expanding  —  about  pf,  the  average  change  in  volume  Aa  at  each 
level,  per  unit  mass,  is: 


.  ^  -IT 

Aa  =  — 

4 


[<fc 

ap 


(Pf) 


+  —  (^) 
dp  ap' 


TT 


da 

dp 


(Pf) 


dp 


& 

dp 


•  it] 

Pf 


-it  d  ,da. 

2  dp  dp 


(1.5) 


Integrating  Aa  over  pressure  as  an  estimate  of  APE  for  this 
simple  case: 


APE (Pf ) 


tt2  da 
2  dp 


ir  da 


(1.6) 


The  right  hand  side  of  (1.6)  is  greater  than  or  equal  to  zero  for  a 
stably  stratified  fluid.  It  will  be  recognized  as  the  Boussinesq 


W 1  tiftf  !!'!!J!i.;'S'M 


24 


approximation  to  APE  per  unit  mass,  and  represents  the  contribution  to 
the  gravitational  potential  energy  from  local  changes  of  the  mass 


field  relative  to  pressure  surfaces  with  no  other  change  resulting 
in  the  mass  field  globally. 

The  small  changes  in  internal  energy  mentioned  earlier  occur 
if  there  are  gradients  of  compressibility  as  a  function  of  vertical 
or  horizontal  position,  since  the  adiabatic  displacement  at  two 
different  positions  in  space  will  cause  the  two  parcels  involved  to 
compress  (or  expand)  differentially,  resulting  in  an  additional  change 
in  the  horizontally  averaged  specific  volume  (Fig.  1.1c).  The  sign 
of  that  average  change  in  volume  depends  upon  the  spatial  gradients 
of  compressibility.  Compressibility  is  primarily  a  function  of 
temperature  (with  colder  water  more  compressible) ,  so  that  these  volume 
changes  may  be  thought  of  as  conversions  between  internal  and  potential 


energy.  That  conversion  process  is  entirely  separate  from  the  local 
changes  of  mass  relative  to  a  pressure  (or  geopotential)  surface.  The 
change  in  volume  due  to  that  conversion  causes  global  changes  in  the  mass 
field.  The  change  in  volume  due  to  compressibility  for  the  upward 
displacement  is,  approximately: 


Act  =  ~K1  (V  Si'  P  +  7T)  ’  * 


and  for  the  downward  displacement: 


AaD  =  k2  (@2,  S2,  p  -  ir)  •  it 


out  n 

with  K  the  adiabatic  compressibility  s'  0  potential  temper- 
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ature  and  S  salinity  (see  Pig.  1.1c).  The  contribution  to  the 
gravitational  potential  energy  per  unit  mass  is 


A  APE 


-/ 


(K.  - 


7T  dp' 


(1.7) 


due  to  horizontal  gradients  of  compressibility.  (A  corresponding  term 
for  vertical  gradients  in  the  general  case  is  derived  as  part  of  the 
rigorous  derivation  given  later  in  this  section.) 

If  the  displaced  columns  are  now  moved  back  to  their  original 

positions  (a  distance  IT) ,  which  corresponds  to  leveling  the  observed 

field,  the  mass  imbalance  caused  by  conversion  of  internal  to  potential 

energy  becomes  obvious  (Fig.  l.ld).  Although  the  figure  illustrates 

only  one  level,  there  are  corresponding  imbalances  at  all  levels.  If 

the  specific  volume  pictured  in  Pig.  l.ld  is  taken  to  be  the  new  'observed' 

field,  af  is  no  longer  the  adiabatic  minimum  at  pf.  Therefore,  in 

order  to  estimate  how  much  of  the  energy  given  by  (1.7)  is  actually 

available  for  conversion  to  kinetic  energy,  the  field  in  Fig.  l.ld  must 

be  leveled,  and  deviations  from  the  new  reference  a,  calculated 

f 

(Fig.  l.le).  The  amount  actually  available  is  a  small  portion  of  that 
calculated  in  (1.7)  since  most  of  the  energy  must  go  into  changing  the 
reference  level  The  energy  given  in  (1.7)  is  itself  a  small  cor¬ 

rection  to  the  Boussinesq  APE  (1.6),  or  order  <  15%  for  most  oceanic 
applications. 

If, rathef  than  considering  gravitational  potential  energy  alone, 
the  sum  of  potential  plus  internal  energies  (enthalpy)  is  considered, 
the  conversions  between  internal  and  potential  energy  will  not  appear 
explicitly.  (Lorenz,  1955  and  Reid,  et  al . ,  1980  both  define  APE  as  the 
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difference  in  enthalpy  between  the  observed  state  and  the  adiabatic 
minimum.)  Since  the  conversion  terms  are  effectively  very  small,  using 
enthalpy  rather  than  gravitational  potential  energy  will  result  in  a 
good  approximation  to  the  true  APE.  In  order  to  use  the  difference 
in  enthalpy,  one  must  either  determine  from  observations  the  reference 
enthalpy  field  or  else  use  an  expansion  about  the  displacement  TT  to 
evaluate  the  change  in  enthalpy  between  the  observed  and  reference 
fields.  The  first  approach  presents  some  practical  difficulties,  but 
could,  in  principle,  be  used.  The  second  approach  is  satisfactory 
unless  the  displacements  are  large,  in  which  case  higher  order  terms 
in  the  expansion  must  be  included.  In  that  case  a  better  approxi¬ 
mation  to  the  true  APE  is  the  gravitational  potential  energy,  which 
can  be  evaluated  exactly,  with  the  small  conversion  terms  like  (1.7), 
whose  errors  for  large  displacement  are  also  small,  subtracted  out. 

In  order  to  make  the  rigorous  thermodynamic  argument  which  cor¬ 
responds  to  the  example  just  given,  a  few  definitions  must  be  made.  In 
the  interests  of  brevity  the  different  estimates  of  APE  are  given 
abbreviations  as  well: 

a.  Available  Gravitational  Potential  Energy  (GPE) 

rr  - 

GPE  (per  unit  mass)  =  /  (Ct.  -  a^)  dp  (1.8) 

P 

(with  and  af  the  initial  and  reference  specific  volumes,  respec¬ 

tively,  p  the  reference  pressure,  and  the  overbar  indicating  hori¬ 


zontal  average) . 
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b.  Total  Available  Gravitational  Potential  Energy  (TGPE) : 


TGPE  (per  unit  area) 


"/  ?/  (0ti  “  af}  dp"  dp> 

P  P" 


c.  The  Boussinesq  Approximation  to  APE  (APE_) : 

B 


(1.9) 


APE  (per  unit  mass) 
B 


q-  a*  .IT2  =  ^N2  C2 

2  p  2  * 


(1.10) 


(with  7T  the  displacement  in  decibars,  a*  the  vertical  gradient  of 

P 

specific  volume,  N  the  buoyancy  frequency  and  £  the  displacement  in 
meters) . 

d.  Available  potential  energy  as  defined  by  Lorenz  (1955)  and 
Reid  et  al.  (1980) (HPE) : 


HPE  (per  unit  area) 


(hi  -  hf)  dp" 


(1.11) 


(Ik  and  hf  the  initial  and  reference  specific  enthalpies,  respectively) . 
The  relationships  between  GPE,  APE_  and  HPE  can  be  demonstrated 

I> 

most  readily  by  considering  perturbation  expansions  in  displacement 
about  the  adiabatically  leveled  reference  state.  For  any  state  variable 
the  connection  between  the  initial  and  reference  states  may  be 
written : 


^(P)  = 


(P  -  ir" 


p-TT 


,3^  „ 

3p  a  dP 


(1.12) 


7f"  =  1T(p  -  TT^) 
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For  =  a,  the  specific  volume,  the  expansion  of  GPE  is  straight¬ 
forward  : 


(0L±  -  af)  dp" 


+  j  a* 

2  p 


—  r  tt2 

(a*ir  +  —3-)  dp" 

r 


(1.15) 


The  first  term  on  the  right  is  the  Boussinesq  APEg  (the  term  at  p^_  is 
evaluated  and  used  as  a  boundary  condition) ;  the  remaining  integral 
term  represents  conversions  of  internal  to  potential  energy.  Physically, 
these  come  about  because  of  the  small  change  in  volume  which'  results 
from  interchanging  parcels  of  water  with  different  compressibilities. 

As  described  earlier,  leveling  moves  colder  (more  compressible)  water 
to  lower  pressure  and  warmer  (less  compressible)  water  to  higher  pressure. 
The  small  change  in  volume  which  results  causes  a  shift  in  the  mass  field 
above  the  location  of  the  leveled  surface.  The  first  half  of  the 
integral  term  represents  conversions  due  to  horizontal  gradients  of  com¬ 
pressibility:  since  daf/dp  is  not  a  function  of  horizontal  position, 

ct*Tr  reduces  to  -KTT  —  the  horizontal  correlation  of  adiabatic  c.,. 

P 

pressibility  and  displacement.  The  second  half  of  the  integral  term 
represents  changes  due  to  vertical  gradients  of  compressibility. 

To  understand  the  relationship  between  TGPE  and  HPE,  an  expansion 
usit  g  \p  =  h  is  suggested: 


+ 


1  _d_ 

2  dp 


h*TT2 

P 


1  _d_ 

2  dp 


e*Tf2 

P 


r  (h)ir2 

_  _ 

2 

2 


+  0  (TT3 ) 

+  0(ir3) 


and 
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with  e  the  specific  internal  energy.  Using  the  following  definition 
of  enthalpy  (Fofonoff,  1962) : 


dh  =  T  dr)  +  yds  +  cxdp  =  d(ap)  +  de 

with  T  absolute  temperature,  T)  the  specific  entropy,  and  y  the 
chemical  potential  of  salinity  S)  the  following  expressions  may  be 
written  in  terms  of  the  specific  volume  Cl: 


rK(h)  "  dp  “  (3p). 


rK(e)  =  a*  -  PrK(a) . 


Then  _ _ 

Ah  -  Ae  =  -P  a*ir  +  ^(p  a*n2) 


pr  (a)  it2 

— ~ - +  0(TT3) 


which,  except  for  a  constant  of  integration,  is  equivalent  to  GPE. 

Thus  HPE  is  the  sum  of  TGPE  and  any  changes  of  internal  energy,  to 

second  order  in  ft.  Since  it  was  shown  earlier  that  GPE  is  the  sum  of 

APE  and  changes  of  potential  energy  due  to  conversions  from  internal 
H 

energy,  it  follows  that  g(dHPE/dp)  and  APE  are  equivalent  to  second 

order  in  IT,  provided  that  APE  is  calculated  using  the  adiabatic  dis- 

B 

placements  from  the  reference  field,  rather  than  the  mean  field. 
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1.4  Procedure  for  calculating  the  adiabatically  leveled  steric  field 

There  are  two  stages  in  the  computational  procedure  for  determin¬ 
ing  the  adiabatically  leveled  reference  field.  In  the  first  stage,  for 
each  CTD  profile,  pressure  and  potential  temperature  0  (Fofonoff,  1977) 


rPf 


9(p»pf) 


=  T(p)  + 


(1.16) 


are  fit  to  nth  order  polynomials  against  potential  specific  volume 
(steric)  anomaly: 


<Hp,Pf)  =  a  (pf  ,0  (p,pf)  ,S  (p))  -  a  (p,0, 35)  .  (1.17) 


The  regressions  are  performed  over  an  interval  of  pressure  Ap  about  pf 
(so  that  0  and  6  are  referred  adiabatically  to  pf  rather  than  to 
atmospheric  pressure)  for  a  number  of  levels.  The  intervals  Ap  may 
overlap. 

in  the  second  stage,  the  coefficients  of  the  pressure  polynomials 
at  each  level  are  averaged  about  a  common  origin  6  : 


6 

m 


1 

M 


(1.18) 


6  the  average  6  over  Ap  for  station  k.  If  the  pressure  regressions 
for  individual  stations  at  a  given  level  are: 


A 


l  a.  (6-6  ) 

u  kn  k 


n=o 


n 


(1.19) 


•|C  I 

MftL  :  I 


then  the  average  coefficients  (about  6  )  define  the  horizontal  average 

m 

of  Pk(6f): 


,  M  ,  M  N 

k=l  k=l  n=o 


(1.20) 


r  (n+r) ! *  4r 

a'  =  )  a,  ,  — —“Co  “0, ) 

kn  L  k,n+r  ntrl  m  k 
r=o 


(1.21) 


The  constraint  of  mass  conservation  during  leveling  requires  the  initial 
mass  p^/g  and  the  final  mass  p^/g  above  the  surface  corresponding  to  6f 
to  be  the  same,  or 


M  N 


>f  ^«£»  -5  X  l  “kn'V6/  • 


(1.22) 


k=l  n=o 


This  equation  may  be  inverted  to  obtain  6^.  Once  6^  is  determined,  the 
displacements  tt,  and  the  initial  and  final  potential  temperatures  0^k 
and  0£  are  given  by 


tt  =  p  -  p  =  l  a  (6  -6  )  -  p 

k  xk  f  L  kn  f  k  f 

n=o 


(1.23) 


ei*=  l  fWW" 


8f=  l  Bk»<W” 

n=o 


(1.24) 
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with  3^^  the  potential  temperature  regression  coefficients  for  a  given 
station.  Note  that  salinity  is  determined  implicitly. 

In  order  to  illustrate  in  a  simple  way  what  the  computational  pro¬ 
cedure  actually  does,  consider  the  case  N  =  1,  i.e.,  specific  volume  a 
linear  function  of  depth,  over  an  interval  Ap  centered  about  pf.  If 
several  stations  are  included,  all  with  6  varying  linearly  with  depth, 
but  with  different  slopes  and  different  average  values  of  6  (referred 
adiabatically  to  p^) ,  then  the  coefficients  c^0  (which  are  equal  to  p^) 
are  the  same  for  all  k,  but  and  <5^  differ  from  station  to  station. 
The  corrected  coefficients  for  averaging  are  (eq.  1.21): 


a*  =  a..  +  a.,(6  -6, ) 
kO  kO  itl  m  k 


akl  akl  ' 


Averaging  over  all  k,  subject  to  the  mass  conservation  constraint  (1.22) 


p£  ’  “o  +  “kl'VV  +  \l(S£-Sm> 


6=6  - 
f  m 


kl  m  k' 


(If  a  is  the  same  for  all  k,  then 

KX 


6=6  . 
f  m 
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For  that  simple  case,  6^  is  just  the  average  of  6  over  Ap  and  all 
stations . ) 

The  value  of  6  at  p.  in  the  observed  field  (6., )  is  just  6,  ;  the 

f  lk  J  k 

contribution  to  GPE  from  p^  is 


TT  "  ’ST(Pf)  '  6f<pf)  ■  \  '  sm  +  “kA-V^l 


■  “ki'W'V 


If  is  a  constant  for  all  k,  the  contribution  to  GPE  is  identically 
zero  (as  expected  from  the  arguments  forwarded  in  the  discussion  of 
potential  energy) .  Notice  that  the  compressibility  effects  have  been 
included  by  referring  the  steric  anomaly  5  to  pf .  The  adiabatic  dis¬ 
placements  7Tk  are  the  difference  between  p^  (the  pressure  corresponding 
to  6  in  the  reference  field)  and  p  (the  pressure  corresponding  to  6 

X  lK  X 

in  the  observed  field  at  station  k) : 


Pik 


=  a 


kO 


+  a  (6  -5  ) 
kl '  f  k 


so  that 


“  \i(  W  • 


Again,  for  the  simple  case  a  ^  a  constant 


7T  =  a  (6  -6  )  . 
k  1  m  k 


ci(5  *“  i 

The  coefficient  is  {— )  ;  TT^  is  then  the  adiabatic  generalization 

(since  6  is  referred  adiabatically  to  pf)  of  the  usual  Boussinesq  dis- 
★ 

placement  TT^: 


* 


~ a 

da 

dp 


In  the  general  case  of  arbitrary  N,  the  higher  order  polynomial  re¬ 
gressions  resolve  the  vertical  structure  of  6  within  Ap,  and  thereby 
provide  an  accurate  weighting  of  the  specific  volume  information  as  a 
function  of  pressure  for  the  specialized  averaging  and  subsequent  in¬ 
version  which  results  in  the  reference  steric  anomaly  at  a  single 
pressure. 

Finally,  the  dynamic  height  AD^,  now  defined  using  af  rather  than 
a(p,0,35)  (Fofonoff ,  1962),  and  GPE  are  calculated  by  numerical  inte¬ 
gration  of 


AD. 


■{ 


<Wd*’ 


(1.25) 


and 


GPE 


(1.26) 


Discussion  of  the  errors  resulting  from  measurement,  finestructure  and 


numerical  errors  is  deferred  to  Chapter  2. 
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■  Chapter  2 

2.1  Introduction 

The  observation  of  low  frequency  signals  in  the  ocean  requires  a 
combination  of  long  time  series  and  adequate  spatial  coverage.  Over  the 
last  twenty  years  the  technology  for  obtaining  long  current  meter  records 
in  the  deep  ocean  has  progressed  sufficiently  that  there  are  now  available 
2  to  3  year  records  which  are  nearly  continuous.  However,  these  measure¬ 
ments  are  necessarily  limited  in  their  horizontal  coverage,  due  to  the 
cost  of  instrumentation.  An  alternative  approach  to  the  problem  is  the 
use  of  hydrographic  measurements  (temperature,  salinity  and  pressure) 
which  have  better  spatial  coverage  both  vertically  and  horizontally,  but 
are  not  continuous  in  time  at  any  level.  From  repeated  hydrographic 
.  ,  measurements  within  a  given  area  over  a  sufficiently  long  period,  however, 

time  series  of  horizontally  averaged  quantities  can  be  constructed,  and 
the  structure  of  low  frequency  motions  can  be  studied.  For  seasonal  sig¬ 
nals  this  requires  a  minimum  of  four  samples  per  year  taken  over  a  period 
of  at  least  a  year,  and  over  a  fairly  extensive  area  whose  dimensions  are 
determined  by  the  dominant  length  scales  of  the  seasonal  signal.  Anti¬ 
cipating  the  results  of  Chapter  3,  the  minimum  wavelength  is  about  400  km. 
L  The  horizontal  resolution  of  the  data  should  be  less  than  a  quarter  of  the 

%  dominant  wavelength,  or  station  spacing  of  no  more  than  about  100  km  for 

i  the  minimum  wavelength. 

One  set  of  hydrographic  data  which  meets  these  requirements  is  a 
l  -  group  of  approximately  600  stations  occupied  by  French  investigators 

|  (Fruchaud,  1975;  Fruchaud,  et  al. ,  1976a,  1976b)  over  a  period  of  three 
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years  in  the  Bay  of  Biscay  off  the  western  coast  of  Prance.  The  data 
consists  of  eleven  cruises  spaced  roughly  three  months  apart,  of  which 
eight  used  CTD  (conductivity,  temperature,  depth)  instruments,  and  three 
used  STD  (salinity,  temperature,  depth)  instruments.  The  data  were  ob¬ 
tained  through  the  United  States  National  Oceanographic  Data  Center,  in 
the  format  in  which  it  was  issued  by  the  French  Organization  Bureau 
National  de  Donnees  Oceaniques.  The  cruises  are  summarized  in  Table  2.1. 

The  majority  of  stations  were  occupied  in  a  region  bounded  by  2°  to 
12°W  and  43°  to  48°N  (see  Fig.  2.1  and  Table  2.1).  For  the  purpose  of 
comparison,  only  stations  within  those  limits  were  used  in  the  calculation 
of  lateral  fields.  The  bottom  topography  of  that  region  (Fig.  2.1)  is 
predominantly  abyssal  plain,  although  a  sharp  shelf  break  and  steep  con¬ 
tinental  slope  mark  the  eastern,  southern,  and  angled  northern  boundaries. 
The  simplicity  of  the  topography  does  not  extend  much  beyond  the  western 
boundary  of  the  smaller  region  chosen  for  study. 

The  original  version  of  the  data  is  an  uneven  pressure  series  of 
approximately  one  decibar  resolution.  Corrections  for  calibration  errors 
were  made  by  the  acquisition  group  (Fruchaud,  1975).  The  original  CTD 
data  was  converted  to  salinity,  temperature  and  pressure  by  the  author, 
using  the  algorithm  from  Fofonoff,  Hayes  and  Millard  (1974).  The  pressure, 
temperature  and  salinity  data  were  then  pressure  sorted  and  smoothed  into 
even  10  decibar  series  in  preparation  for  the  adiabatic  leveling  calcu¬ 


lation. 
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TABLE  2.1 

Description  of  Phygas  Cruises 


CRUISE 

DAY1 

(from  1  Jan  1972) 

DEPTH2 

INSTRUMENT 

TYPE 

NO.  S 

STATIONS 

1. 

Phygas 

22 

124 

- 

135 

1200 

CTD 

31 

(15) 

2. 

Phygas 

23 

235 

- 

242 

1250 

CTD 

70 

(21) 

3. 

Phygas 

24 

294 

- 

307 

1250 

CTD 

68 

(24) 

4. 

Phygas 

31 

383 

- 

393 

1750 

CTD 

58 

(37) 

5. 

Phygas 

32 

480 

- 

492 

1750 

CTD 

60 

(38) 

6. 

Phygas 

33 

554 

- 

567 

1700 

CTD 

44 

(33) 

7. 

Phygas 

34 

627 

- 

641 

1300 

STD 

76 

(50) 

8. 

Phygas 

41 

780 

- 

794 

1750 

STD 

54 

(41) 

9. 

Phygas 

42 

908 

- 

920 

1800 

CTD 

62 

(48) 

10 

.  Phygas 

43 

982 

- 

994 

1750 

STD 

39 

(39) 

11 

.  Phygas 

44 

1076 

- 

1097 

1750 

CTD 

47 

(44) 

Covers  the  period  during  which  stations  (within  the  small  box)  which 
extend  deeper  than  1000  db  were  taken. 

Maximum  depth  common  to  at  least  10  stations,  after  regressions, 
for  the  restricted  region  (Maximum  depth  before  regressions:  add 
250  dbar.) 


3 


Total  number  of  stations.  In  parenthses:  number  of  stations  with¬ 
in  the  restricted  region  (2°-12°W,  43°-48°N)  deeper  than  250  m. 


General  location  of  Bay  of  Biscay  stations,  showing 
bottom  topography.  The  larae  box  encloses  all  stations 
taken  during  the  tree  year  period;  the  smaller  box 
outlines  area  of  greatest  concentration  of  stations. 

The  center  of  the  smaller  box,  referred  to  in  the  text 
as  the  origin,  is  marked  with  a  cross. 
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2.2  Application  of  the  adiabatic  leveling  technique 

The  variable  parameters  used  in  the  initial  regressions  p(6), 

0  (6)  are  summarized  in  Table  2.3,  with  the  values  which  were  used  in  all 
eleven  cruises.  This  initial  calculation  was  performed  for  all  stations 
in  all  cruises.  The  horizontal  averaging  and  determination  of  the  refer¬ 
ence  specific  volume  field,  the  displacements,  and  the  potential  temper¬ 
ature  (hence  implicitly  salinity)  on  the  reference  and  observed  specific 
volume  surfaces  was  performed  separately  for  the  subset  of  stations  with¬ 
in  2°-12°W,  43°-48°N,  for  each  cruise;  the  results  which  follow  are  only 
those  from  this  restricted  region.  The  number  of  stations  in  the  restric¬ 
ted  region  for  each  cruise  is  given  in  Table  2.1. 

There  are  several  sources  of  uncertainty  in  the  adiabatic  leveling 
calculation,  resulting  from  measurement,  finestructure,  and  numerical 
errors.  Methods  for  calculating  the  corresponding  errors  in  derived 
quantities  are  described  below.  The  errors  themselves  will  be  found 
plotted  in  various  figures  and  listed  in  tables  as  noted.  For  a  complete 
discussion  of  the  relative  errors  in  the  different  estimates  of  APE, 
the  reader  is  referred  to  Bray  and  Fofonoff  (1980) .  Numerical  errors 
are  only  important  in  calculations  involving  numerical  integration 
(GPE,  TGPE,  HPE,  dynamic  height) ,  and  result  from  inadequate  resolution 
of  vertical  structure.  A  method  for  determining  the  correct  vertical 
resolution  of  the  adiabatic  leveling  technique  in  order  that  integration 
errors  are  made  smaller  than  non-reducible  errors  is  described  in  detail 
also  in  Bray  and  Fofonoff  (1980) . 

This  discussion  concentrates  on  the  effects  of  measurement  and 
finestructure  errors  on  derived  quantities  such  as  vertical  displacement, 
GPE,  and  maps  of  salinity  or  temperature  on  leveled  steric  surfaces.  By 


TABLE  2.3 


Initial  Calculation  -  Regression  Parameters 


pf 

step 

Ap 

0,200 

(50) 

250 

250,350 

(50) 

400 

400,450 

(50) 

500 

500,1100 

(50) 

450 

1150,2000 

Knn 
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measurement  error  is  meant  the  uncertainty  in  the  measured  values  of 
pressure,  temperature,  and  salinity  (or  conductivity)  which  results 
from  uncorrected  calibration  errors  (quantization  noise  is  ignored)  and 
which  will  be  taken  to  be  randomly  distributed  from  one  station  to  the 
next,  but  constant  for  a  given  station;  random  finestructure  errors  are 
the  fluctuations  caused  by  real  variability  of  the  field  which  occurs  on 
smaller  scales  than  those  over  which  the  calculations  are  performed.  Both 
of  these  errors  will  be  treated  as  random  errors  in  this  discussion;  the 
values  for  the  variance  of  pressure,  temperature,  and  salinity  due  to 
measurement  error  are  taken  from  Fruchaud  (1975)  (see  Table  2.2).  As  es¬ 
timates  of  the  vertical  finestructure  errors,  values  the  variance  of 
the  regression  estimates  over  each  regression  interval  Ap  are  calculated 
from  the  regression  residuals,  following  Fofonoff  and  Bryden  (1975) .  An 
estimate  o  f  the  variance  02  of  a  single  data  point,  for  either  temper¬ 
ature  or  pressure  is: 


Q 


V=1 

with  Q  the  number  of  observations  in  Ap,  p  .  Q  the  observed  pressure 

V  V 

A  A 

and  potential  temperature,  and  p^,  ©v  the  regression  estimates  of  p 
and  0  at  6^  .  Further,  estimates  of  the  variance  of  the  regression 
estimates  p,  0  at  any  value  of  6  are  given  by 
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N+l  N+l 

v<P«)>  -  al  I  I 

i=l  j=l 


(6  -  6  )1+j“2 
it 


v(0  (6) )  =  a 


=  ^2 
6 


N+l  N+l 

II 

i=l  j=l 


-1 


ID 


«$  - 


V 


i+j-2 


(2.2) 


R.  . 

3-D 


■I 


(6..  “ 


V 


i+j-2 


V=1 


Given  these  estimates  for  measurement  and  finestructure  errors  in 
pressure,  temperature  and  salinity,  errors  in  derived  quantities  such  as 
steric  anomaly  and  Brunt-Vaisala  frequency  can  be  calculated  in  the  same 
way  as  in  Gregg  (1979) ,  using  Taylor  series  expansions  in  the  independent 
variables  p,  T  and  S.  The  errors  in  the  vertical  displacements  ir^: 


are  given  by  the  errors  in  p^,  since  pf  is  chosen  and  therefore  has 
no  error  associated  with  it.  (The  errors  in  p^  are  just  the  pressure 
errors.)  Thus  the  errors  in  each  term  of  the  expansion  (1.15)  of  GPE 
can  be  estimated  as  a  function  of  the  displacement  error,  as  described 
in  detail  in  Appendix  A. 

The  reference  specific  volume  surfaces  6f  have  small  random  errors 
because  of  the  heavy  vertical  and  horizontal  averaging  which  is  inherent 
in  the  leveling  technique.  The  polynomial  regressions  of  p  and  0 
against  steric  anomaly  allow  an  accurate  weighting  of  the  vertical 
average,  without  impairing  vertical  resolution.  The  measurement  errors 
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in  6.  are: 
f 


— 

°  r» 


/1£\  2  2  ,  ,9£\2  e2  +  /M\2  p2 
9T  £T  +  V  %  +  (9S  CS 


E  ,  e  ,  e  the  measurement  errors  in  T,  p,  and  S  respectively. 

jp  *3 

(Recall  that  measurement  errors  are  taken  to  be  random  here.)  Random 
errors  in  0f  and  Sf  mapped  on  6f  surfaces  are  also  small,  since 
0  (and  implicitly  S)  are  known  as  functions  of  6;  hence  random  errors 
in  0 and  Sf  are  estimated  by  the  product  of  the  averaged  gradients 
(90/9 9S/9Sf)  and  the  small  random  errors  in  6^.  These  two  types 
of  averaging  (for  and  for  0^*  separately)  combine  to  reduce  sig¬ 

nificantly  the  random  errors  in  0  and  S  mapped  on  reference  steric 
surfaces 

As  discussed  by  Gregg  (1979)  other  less  heavily  averaged  techniques 
for  mapping  0  and  S  on  density  surfaces  can  introduce  large  random 
errors  due  to  f inestructure .  Uncorrected  bias  errors  resulting  from  the 
use  of  more  than  one  CTD  can  also  introduce  large  errors.  The  latter 
are  reduced  by  the  adiabatic  leveling  technique,  since  0  and  S  are 
computed  as  functions  of  6,  rather  than  p,  and  the  averaged  gradients 
of  0  and  S  with  6  are  less  variable  than  the  local  gradients  with  p. 
For  example,  if  the  standard  method  outlined  by  Gregg  is  used  to  map  0 
and  S,  then  the  bias  error  in  potential  temperature  can  be  estimated  as 


A0  =  A  -(||)  ,>5) 

6  8p  local'  9p 


local 


whereas  the  estimated  bias  error  in  0^  is 


V 
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( 


(Ag  the  bias  error  in  specific  volume) .  Gregg  points  out  that  it  is 
primarily  the  local  gradients  of  9,  S,  and  6  with  pressure  which 
cause  the  errors  in  0  and  S  to  be  large  when  mapped  on  density  sur¬ 
faces;  that  problem  does  not  exist  in  the  determination  of  0f  and  S^. 
Furthermore,  the  regressions  of  0  and  p  provide  accurate  estimates 
of  the  vertical  finestructure  errors  in  the  observed  field.  The  deter¬ 
mination  of  horizontal  finestructure  errors  in  0  and  S  mapped  on 
reference  steric  surfaces  requires  more  information  than  is  given  by  the 
adiabatic  leveling  technique  alone.  In  the  next  section  a  method  is 
described  for  determining  those  errors  in  a  fashion  analogous  to  that 
given  earlier  for  vertical  finestructure. 

The  errors  in  derived  quantities  described  above  generally  vary 
with  depth;  typical  values  are  given  in  Table  2.4,  and  more  precise 
estimates  are  plotted  as  error  bars  in  the  appropriate  figures. 


50 

TABLE  2.4 

Typical  error  values  for  derived  variables 


Variable 

Measurement 

Error 

Vertical  Finestructure 
Error 

ir 

±5  db 

±5  db 

6f 

±0.5  x  10  5  cm3  •  gm  1 

±0.2  x  lcT5 

cm3 

•  gm  1 

S(6f) 

±.010  ppt 

±.010  ppt 

0(6f) 

±.010°C 

±.010°C 

Ad 

±.l  dyn  cm 

±.05  dyn  cm 

GPE 

±15.  x  lo"1*  j  •  kg"1 

±10.  x  lo"4 

J  • 

kg"1 

TGPE 

±0.02  x  lo"4  3  •  m”2 

±.01  x  lo'4 

J  • 

m~2 
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2.3  Mean  fields  and  non-seasonal  variability 

As  a  background  to  the  discussion  of  seasonal  variability,  a 
description  of  the  time  mean  fields  in  the  region  2°  to  12°W,  43°  to 
48°N  is  presented  here.  The  discussion  is  divided  into  three  parts: 
the  first  deals  with  variables  which  provide  a  basic  description  of  the 
physical  system  (reference  steric  field,  buoyancy  frequency,  potential 
temperature-salinity  relationship) ;  the  second  with  variables  related 
to  the  energetics  of  the  system  (vertical  displacements,  dynamic  height, 
APE) ;  and  the  third  with  variables  which  may  be  used  as  tracers  (salinity 
anomaly,  temperature  variance) . 

The  vertical  structure  of  the  reference  steric  field  changes 

very  little  with  time.  A  typical  plot  of  6^  vs.  pressure  is  shown  in 
Fig.  2.2.  There  is  a  deep  nearly  mixed  layer  extending  to  about  400  m, 
with  a  weak  main  pycnocline  below,  and  very  little  change  in  below 

the  pycnocline.  The  average  buoyancy  frequency  (N)  profile  in  Fig.  2.3, 
which  was  constructed  by  averaging  N  horizontally  at  each  pressure 
p^  for  each  cruise,  and  then  averaging  cruises,  shows  the  same  structure 
in  a  different  perspective.  There  is  a  seasonal  thermocline  which  is 
not  plotted.  (The  adiabatic  leveling  technique  is  not  valid  shallower 
than  about  300  m,  since  density  surfaces  are  not  continuous  horizontally, 
but  may  intersect  the  surface  of  the  ocean.  All  levels  shallower  than 
300  m  have  thus  been  excluded  from  consideration  here.)  The  average  6 
vs.  S  relationship  is  plotted  in  Fig.  2.4,  along  with  the  Worthington 
and  Metcalf  (1951)  0  vs.  S  standard  curve  for  the  western  North  Atlan¬ 

tic  for  comparison.  The  dominance  of  Med  Water  influence  below  the  main 


thermocline  is  clear. 
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Figure  2.3  Buoyancy  frequency  N,  averaged  over  all  stations,  with 
the  range  of  values  from  individual  cruise  averages 
given  by  horizontal  lines  at  each  depth. 


Average  9  vs  S  for  all  cruises  (solid  line).  Dotted 
line  is  a  cubic  spline  fit  to  6  vs  S  from  Worthington 
and  Metcalf  (1961)  and  Iselin  (1939)  for  the  western 


North  Atlantic. 
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The  root  mean  square  (rms)  vertical  displacements  for  each  pressure 
pf  and  cruise  are  shown  in  Fig.  2.5.  Neither  the  average  stratification 
nor  the  rms  displacements  bears  a  strong  resemblance  to  the  typical  plot 
of  GPE  shown  in  Fig.  2.6,  which  has  a  strong  second  maximum  below  the 
main  thermocline.  However,  a  plot  of  GPE  with  vertical  structure  similar 
to  Fig.  2.6  can  be  constructed  by  multiplying  the  average  N2  and  C2 
(displacement  in  meters)  at  each  level  as  shown  in  Fig.  2.7.  This  is 
a  Boussinesq  type  of  estimate,  except  that  correlations  between  t; 
and  N  are  neglected,  so  that  the  amplitude  which  results  differs  from 
that  of  a  typical  plot  of  GPE.  It  is  the  combination  of  a  slowly 
decreasing  stratification  and  increasing  rms  displacement  which  causes 
the  second  maximum  in  GPE  below  the  thermocline.  For  a  comparison  with 
western  North  Atlantic  GPE  estimates  the  reader  is  referred  to  Bray 
and  Fofonoff  (1980)  who  discuss  the  application  of  the  adiabatic 
leveling  technique  to  the  Mid-Ocean  Dynamics  Experiment  (MODE) .  The 
second  maximum  in  GPE  does  not  appear  in  any  of  the  MODE  GPE  profiles. 
Estimates  of  eddy  kinetic  energy  (KE)  from  a  current  meter  mooring 
located  at  47°N,  10°W  as  part  of  the  NEAD  observational  program  are  also 
plotted  in  Fig.  2.6  (Colin  de  Vediere,  personal  communication).  In  the 
thermocline  GPE  exceeds  KE  by  a  factor  of  3;  because  of  the  second  maximum 
in  GPE  below  the  thermocline,  GPE  at  1500  db  exceeds  KE  there  by  a  factor 
of  8,  suggesting  that  there  is  storage  of  potential  energy  in  the  density 
field.  The  vertical  integral  of  GPE  (TGPE)  between  300  and  1100  db  has 
a  mean  value  over  all  cruises  of  0.500  x  10  in  2  (Table  2.5). 


Measurement  and  finestructure  errors  in  TGPE  are  also  given  in  Table  2.5. 
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Figure  2. 


Rms  vertical  displacements  for  all  cruises 


Figure  2.6  APE  for  Phygas  42.  Open  circles  are  GPE,  squares  are 
APE  ,  x  and  +  are  the  contributions  from  horizontal 
and  vertical  gradients  of  compressibility  respectively, 
and  triangles  are  the  sum  of  the  Taylor  expansion  (1.15). 
Error  bars  with  heavy  end  lines  represent  measurement 
errors,  those  with  lighter  lines  represent  finestructure 
errors.  Isolated  solid  circles  at  600,  1000,  and  1500  db 
are  estimates  of  eddy  kinetic  energy  (tides  and  inertial 
motions  removed)  from  a  current  meter  mooring  at  47 °N,  10°W 


designated  NEAD7. 


TABLE  2.5 


Cruise 


Potential  Energy  Per 

Unit  Area 

TGPE 

Ap.e. 

Errors 

Finestructure 

(%) 

Measure 

22 

lO-1*  J  •  m~2 

1.043 

I-. 5 

1 

o 

1 — l 

. 

•  m 

14.0 

10.4 

23 

0.465 

-3.21 

2.2 

4.7 

24 

0.609 

-8.12 

1.3 

3.9 

31 

0.416 

6.10 

2.2 

4.8 

32 

0.485 

5.94 

3.1 

4.5 

33 

0.384 

-7.43 

3.1 

5.2 

34 

0.525 

3.61 

1.5 

4.0 

41 

0.342 

0.48 

5.2 

5.2 

42 

0.409 

1.22 

1.5 

4.9 

43 

0.559 

0.36 

1.6 

4.1 

44 

0.261 

-2.58 

2.3 

6.9 

Mean  (std  dev):  ,500(.206)  -0.36(4.99) 

300  _P 

,  f 

1  AP 


;  /“  /'  [6^  (n+1)  -  6^(n)]  dp'dp,  n  the  sequential 


1100  1100 


number  of  the  cruise 
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Dynamic  height  at  500  db  referred  to  1000  db  calculated  according 
to  (1.25)  is  plotted  for  each  cruise  in  Figs.  2.8  to  2.10.  The  plots 
are  objective  maps  with  an  imposed  correlation  scale  of  100  km,  as 
indicated  by  correlation  function  calculations.  The  signal  in  dynamic 
height  is  weak  (a  maximum  of  ±  3  dyn  cm)  and  generally  correlated  with 
the  dispacements  it  between  500  and  1000  db  (also  plotted  in  the 
dynamic  height  maps) .  The  dynamic  height  field  is  not  steady  but  shows 
evidence  of  slow  westward  propagation  of  features  at  2  to  5  cm  sec  1 
(Fig.  2.11). 

Because  the  changes  in  salinity  over  the  field  are  small  compared 
to  the  salinity  itself,  and  features  of  interest  may  be  obscured, 
salinity  anomaly  (defined  as  the  difference  between  the  observed  salinity 
and  the  salinity  which  a  parcel  of  water  at  the  observed  potential  tem¬ 
perature  would  have  if  it  adhered  to  the  Worthington-Metcalf  and  Iselin 
0— S  curves)  was  examined.  Arrai  (personal  communication)  fitted  a  cubic 
spline  to  the  combined  Worthington  and  Metcalf  (1961)  0-S  curve  for 

the  western  North  Atlantic  and  a  modified  version  of  Iselin' s  (1939) 

T-S  curve  for  the  intermediate  and  warmer  waters,  in  which  temperature 
was  converted  to  potential  temperature  using  standard  temperature-pressure 
correlations  from  hydrographic  data.  The  coefficients  for  this  cubic 
spline  fit,  which  were  modified  slightly  in  the  near  surface  waters  by 
Maillard  (personal  communication)  are  found  in  Appendix  B.  The  curve  is 
plotted  in  Fig.  2.4. 

The  average  over  all  cruises  of  the  horizontally  averaged  salinity 
anomaly  on  reference  steric  surfaces  is  plotted  in  Fig.  2.12.  There  is 
some  change  in  the  average  structure  between  cruises;  however,  the  hori¬ 
zontal  average  for  a  given  cruise  is  biased  by  the  location  of  the  stations 


Figures  2. 
to  2.10 


8  Objective  naps  of  dynamic  height  at  500  relative  to  1000  db 
for  all  cruises  (contour  units  are  1  dyn.  cm.).  The 
horizontal  coordinates  are  kilometers  from  the  origin 
position,  chosen  to  be  45.5°N,  7°W.  .Superimposed 
upon  the  dynamic  height  contours  are  the  displacements 
it  as  a  function  of  depth  for  each  station.  The  symbol 
identifies  the  cruise  and  is  located  at  the  geographical 
position  of  the  station.  It  also  marks  the  1000  db  depth 
on  the  pressure  axis.  The  scales  for  pressure  and  displace¬ 
ment  are  given  in  the  inset  on  each  figure.  The  cruises 
are  divided  up  by  year:  Fig.  2.8  is  1972  (Phygas  22,  23,24), 
Fig  2.9  is  1973  (Phygas  31,  32,  33,  34)  and  Fig.  2.10  is 


1974  (Phygas  41,  42,  43,  44). 


Figure  2-11  Contour  plot  of  dynamic  height  as  a  function  of  time 

(vertical  axis)  along  a  line  200  km  south  of  the  origin 
(45.5°N,  7°W) .  Contour  interval  is  1  dyn.cm. 
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since  there  are  strong  mean  horizontal  gradients,  as  can  be  seen  in  the 
objective  maps  of  salinity  anomaly  on  the  reference  surface  at  1000  dbar 
(Figs.  2.13-2.15). 

The  mean  horizontal  structure  of  the  salinity  anomaly  field  was 
determined  using  a  least  squares  linear  regression  in  x  (east)  and  y 
(north)  distance  from  the  origin  at  45.5°N,  7°W  (see  Fig.  2.1): 

M  N 

5  -  Z  I  Sij  lx  -  xc>i  -  V3  (2-3) 

j=0  i=0 

with  x  and  y  the  origin  coordinates.  These  regressions  were 
o  o 

performed  over  steric  surfaces  corresponding  to  400,  600,  800,  1000  and 
1200  db  for  all  cruises,  and  also  over  1400,  1600,  1700  db  for  cruises 
with  sufficiently  deep  data.  The  optimum  number  of  terms  in  the  regres¬ 
sion  was  determined  by  examining  the  statistical  confidence  of  the  ratio 
of  each  coefficient  As. .  to  its  standard  deviation.  Four  terms  were 

identified  at  95%  confidence  for  most  levels  and  cruises:  AS  ,  the 

oo 

average  value,  AS^,  the  zonal  (x)  gradient:  ASQ1,  the  meridional  (y) 
gradient;  and  AS^,  the  zonal  curvature  or  second  derivative.  The 
mean  zonal  and  meridional  gradients  (using  all  stations)  are  listed  in 
Table  2.6.  Note  that  the  gradients  are  evaluated  at  the  origin. 

Horizontal  finestructure  errors  in  salinity  anomaly  were  estimated 
by  examining  the  variance  of  the  regression  estimate  of  As  as  was  done 
in  section  2  to  determine  the  vertical  finestructure  errors  in  0  and  p. 
The  estimates  of  error  in  As  are  plotted  in  Fig.  2.17.  Similar  regres¬ 
sions  were  done  using  the  displacement  V  and  the  vertical  finestructure 

error  a_.  A  mean  horizontal  gradient  in  TT  would  indicate  a  mean  slope 
0 

in  the  isopycnals;  however,  no  discernable  (i.e.,  95%  confidence)  grad- 


Figures  2.13  Objective  maps  of  salinity  anomaly  at  1000  db  for  all 


to  2.15  cruises.  Contour  internals  are  .050  ppt;  imposed 

correlation  scale  is  100  km. 
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TABLE  2.6 


Gradients  of  Salinity  of  Anomaly  at  Origin 


Mean1 

Mean2 

Number  of 

S  (s  ) 

x  10 

S  (S  ) 
y  01 

realizations3 

Depth 

10s  ppt  •  km  1 

10s  ppt  •  km  1 

S 

X 

S 

y 

400 

0.00  (.03) 

-0.'05  (.01) 

11  (5) 

(11) 

600 

-0.02  (.04) 

-0.09  (.02) 

11  (3) 

(ID 

800 

-0.05  (.16) 

-0.24  (.05) 

11  (3) 

(ID 

1000 

-0.11  (.16) 

-0.34  (.08) 

11  (5) 

(ID 

1200 

-0.12  (.21) 

-0.42  (.13) 

11  (7) 

(11) 

1400 

-0.10  (.11) 

-0.40  (.10) 

7  (1) 

(7) 

1600 

+0.06  (.16) 

-0.28  (.03) 

7  (1) 

(7) 

1700 

+0.08  (.13) 

-0.23  (.04) 

6  (1) 

(6) 

1  Mean 

over  all  cruises  of 

value  of  at 

origin.  In 

parentheses 

standard  deviation  over  all  cruises. 

2  Mean  over  all  cruises  of  value  of  5^  at  origin.  In  parentheses: 
standard  deviation  over  all  cruises. 

3  In  parentheses:  —  number  of  realizations  for  which  is  dif¬ 

ferent  from  zero  at  the  95%  confidence  level;  S  —  same  as  S  , 

y  x 

except  for  . 
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ients  in  it  were  found.  The  finestructure  error  in  potential  temperature 
is  an  indication  of  the  strength  of  vertical  mixing:  in  general  the 
greater  the  finestructure,  the  less  well  mixed  is  the  fluid.  There  are 
strong  zonal  gradients  of  temperature  finestructure  in  these  data,  with 
the  strongest  mixing  occurring  near  the  boundary  (see  Table  2.7).  The 
zonal  gradient  increases  with  depth  between  400  and  1200  db  which  may 
simply  reflect  the  strong  variability  offshore  associated  with  the  Med 
Water  which  has  a  maximum  at  roughly  1200  db.  The  meridional  gradient  of 
temperature  finestructure  is  only  occasionally  discernable  at  95%  con¬ 
fidence  . 
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TABLE  2.7 

Zonal  Gradient  of  Temperature  Finestructure 


Pressure 

Mean  zonal 

gradient 

(10~5  °C  km"1) 

Std.  dev. 

Number  of 

realizations 

400 

3.9 

2.7 

11 

CO 

v/ 

600 

5.3 

1.7 

11 

(9) 

800 

9.6 

3.7 

10 

(8) 

1000 

11.2 

3.0 

10 

(10) 

1200 

13.6 

3.9 

10 

(10) 

1400 

12.0 

6.0 

7 

(4) 

In  parentheses:  number  of  realizations  significantly  different 
from  zero  at  greater  than  95%  confidence.' 


l 


83 


2.4  Seasonal  variability 

There  .  re  two  types  of  seasonal  signal  below  the  main  thermocline 
in  the  Bay  of  Biscay  data:  an  oscillation  of  the  reference  s\ jric  field 
which  necessarily  occurs  on  horizontal  scales  of  at  least  1000  km  (the 
horizontal  extent  of  the  data),  since  6f  is  independent  of  horizontal 
position;  and  seasonal  variations  in  the  horizontal  structure  of  salinity 
anomaly,  which  occurs  on  smaller  horizontal  s cales.  Because  there  are 
only  eleven  data  points  in  the  time  series,  a  simplified  harmonic  analy¬ 
sis  is  used  to  identify  the  frequency,  jr.pl itude  and  phase  of  the  signals. 
Phase  and  amplitude  at  a  given  frequency  to  are  calculated  by  a  least 
squares  regression  of  the  form 

^-(t)  =  a  cos  'wt  +  <j>) 

=  a^  cos  „t  +  a2  sin  wt  (2.4) 

with  a^  =  a  cos  <f> 

a?  -  -a  sin  <j> 

^  any  variable,  <}>  the  phase,  and  cos  tot,  sin  tot  the  independent 
variables.  By  stepping  through  a  range  of  frequency  and  comparing  the 
statistical  confidence  of  the  coefficients  a^  and  the  best  fit 

of  frequency  can  be  established  .  The  ratio  of  each  coefficient  to  its 
standard  de\  '..jJ'ion  for  frequencies  corresponding  to  periods  of  8  to  16 
months  are  plotted  in  Fig.  2.16  for  salinity  anomaly  (ASoq)  at  1000 
The  only  period  for  which  both  coefficients  have  greater  than  95%  con¬ 


fidence  is  12  months. 


Figure  2.16  Ratios  of  the  coefficients  of  the  harmonic  analysis  (2.4) 
to  their  standard  deviations  for  a  range  of  frequency. 
The  values  at  frequency  1  cpy  are  circled. 
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There  are  seasonal  signals  also  at  800  and  1200  db  (Pig.  2.17); 
there  are  too  few  points  below  1200  db  to  identify  a  seasonal  signal 
with  95%  confidence.  The  signal  double  amplitude  (peak  to  trough)  is 
40  ppm  at  800  db  and  50  ppm  at  1000  and  1200  db.  For  the  two  deeper 
levels  the  phases  calculated  correspond  to  maximum  values  at  130  and  134 
days  from  1  January.  The  difference  in  phase  is  less  than  the  probable 
error  in  phase  determination.  The  signal  at  800  db  is  not  strong  enough 
to  allow  a  determination  of  phase,  so  a  phase  of  130  days  was  used  in 
the  regression  and  only  amplitude  calculated.  None  of  the  higher  order 
terms  As10'  ^Sqi'  ^S20  s^ows  a  clearly  seasonal  signal,  although 
AS20'  the  zonal  curvature  term,  does  change  in  time,  unlike  either 
gradient  term.  The  change  in  AS2Q  cann°t  be  clearly  identified  as 
seasonal.  A  separate  analysis  of  horizontal  structure  was  performed  in 
which  the  curvature  (As.^)  term  was  omitted  from  the  horizontal  regres¬ 
sion.  The  AS  terms  which  result  are  different  from  those  of  the 
oo 

initial  calculation  and  show  no  seasonal  variability.  This  is  interpre¬ 
ted  as  evidence  that  the  seasonal  signal  in  AS  results  from  changes 

oo 

in  the  horizontal  structure  of  the  salinity  anomaly  field. 

The  reference  steric  field  also  changes  seasonally:  at  1000  and 

—  5  3 

1200  db  the  double  amplitude  is  2  *  10  cm  gm  with  phase  set  at 
130  days  from  1  January  (see  Fig.  2.18).  This  represents  a  seasonal 
displacement  of  the  reference  steric  surfaces  of  ±  15  m.  There  is  no 
evidence  from  this  analysis  that  there  is  a  seasonal  signal  in  the  ver¬ 
tical  displacements  Ti;  however  there  are  changes  in  the  vertical  struc¬ 
ture  of  GPE  which  may  be  seasonal  (see  Fig.  2.19).  The  strongest  maxima 
(relative  to  the  thermocline)  in  GPE  below  the  thermocline  occur  in  the 
winter  cruises  (3i  and  41)  and  summer  cruises  (33  and  43) .  The  reference 


Figure  2.17  Salinity  anomaly  at  origin  (45.5°N,  7°W)  ,  as  estimated 
by  horizontal  regression,  vs  time.  Error  bars  are 
horizontal  finestructure  uncertainties  in  the  coefficient 
of  the  zeroth  order  polynomial  terra,  AS>00*  Solid  lines 
are  harmonic  analysis  fits  of  amplitude  and  phase  with 
frequency  1  cpy. 
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Figure  2.18  Reference  specific  volume  anomaly  6^  vs  time.  Error 

bars  are  measurement  errors,  which  are  larger  than  fine- 
structure  errors  inferred  from  differences  in  uptrace 
and  downtrace  calculations.  Solid  lines  are  harmonic 
analysis  fits  of  amplitude;  phase  is  the  same  as  for 
1000  db  ASqo  (Fig.  2.17)  and  frequency  is  1  cpy. 


REFERENCE  SPECIFIC  VOLUME  ANOMALY  8f  vs. TIME 
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\  i 


level  (1100  db)  GPE  (Fig.  2.19,  inset)  also  appears  to  vary  with  seasonal 

periodicity  but  that  signal  cannot  be  identified  with  95%  confidence. 

Of  primary  interest  are  the  horizontal  scales  of  the  different 

seasonal  signals,  in  particular,  that  of  the  salinity  anomaly  signal. 

A  rough  estimate  of  that  scale  can  be  made  by  comparing  the  time  changes 

of  the  zero  order  term  (AS  )  and  the  zonal  curvature  term  (As  ) : 

oo  20 

_3_ 

l2  ^  3t  Aoo  ^  25  x  IQ  3  ppt 

3  AC  5  x  10  7  ppt/km2 

3t  A  20 

or 

a, 

L  =  200  km. 

However,  if  the  signal  is  actually a-propagating  wave,  it  cannot  be  iden¬ 
tified  properly  using  this  technique.  Given  a  zonal  wavelength,  the 
data  can  be  fitted  to  a  propagating  wave  model,  with  amplitude  and  phase 
determined  by  least  squares  techniques.  The  new  phase  would  then  contain 
only  time  information,  not  mixed  time  and  space  information  as  it  does 
in  this  analysis, 

In  the  next  chapter  a  model  is  constructed  in  which  the  wave  lengths 
of  the  dominant  free  propagating  waves  are  calculated,  and  the  resulting 
change  in  salinit.  inomaly  estimated.  The  model  predicts  seasonal  vari¬ 
ability  in  the  displacement  field  as  well.  Using  the  model  prejudice 
for  horizontal  scale,  the  data  are  refitted  to  traveling  plane  waves  with 
seasonal  period  and  good  agreement  is  found  between  the  model  and  the 
data  in  the  salinity  anomaly  signal  and  also  in  the  displacement  field, 
which  is  found  to  have  a  seasonal  signal  when  the  appropriate  analysis 


is  applied. 


Figure  2.19  GPE  for  all  cruises,  with  the  reference  value  (1100  db) 
subtracted,  for  comparison  of  vertical  structure.  Inset 
has  APE  (1100  db)  vs  time. 


relative 
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2.5  Discussion 

From  this  analysis  of  the  Bay  of  Biscay  observations  the  following 
conclusions  may  be  drawn: 

1)  There  is  a  well  defined  mean  field  in  salinity  anomaly  in  the 
Bay  of  Biscay  between  depths  of  400  and  1800  m.  The  horizontal  structure 
is  dominated  by  a  strong  meridional  gradient  with  a  maximum  value  of 
-.42  x  lo  3  ppt  km  1  at  1200  db;  there  is  a  weaker  zonal  gradient  with 
maximum  value  of  -.12  x  10  3  ppt  km  1  at  the  same  depth  and  zonal  cur¬ 
vature  with  a  maximum  value  of  -.5  x  lo  6  ppt  km  2  also  at  1200  db. 

2)  Departures  from  the  mean  horizontal  structure  below  the  thermo- 
cline  are  seasonal  and  result  from  changes  in  the  zonal  structure,  evi¬ 
denced  by  changes  in  the  zonal  curvature  term.  The  resultant  change  in 
salinity  anomaly  at  a  given  location  is  0(±  25  ppm) . 

3)  The  reference  specific  volume  field  oscillates  seasonally 
with  an  amplitude  of  about  1  x  10  5  cm3  gm  1 ,  which  corresponds  to  a 
vertical  displacement  of  the  thermocline  of  ±  15  m.  The  scale  of  this 
motion  is  not  determined  by  the  observations,  but  must  be  larger  than 
the  spatial  extent  of  the  data,  since  the  reference  field  is  independent 
of  horizontal  position. 

4)  The  vertical  structure  of  APE  per  unit  mass  consists  of  two 
maxima:  one  at  the  thermocline,  and  one  of  the  same  approximate  magnitude 
at  about  1400  db.  The  thermocline  APE  is  of  the  order  of  7o  x  10  k  J  kg 
The  amplitude  and  the  vertical  structure  of  APE  change  between  realiza¬ 
tions.  That  change  is  not  demonstrably  seasonal,  although  there  is  visual 
similarity  between  cruises  which  occurred  at  the  same  time  of  the  year. 

The  unresolved  questions  raised  by  the  analysis  of  observations  are: 
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1)  What  are  the  horizontal  scales  of  the  seasonal  variations  in 
salinity  anomaly  and  density? 

2)  What  causes  the  observed  seasonal  variations? 

3)  Should  corresponding  signals  be  expected  in  vertical  displace¬ 
ment  and  hence  in  APE? 

In  Chapter  3  a  theoretical  model  of  oceanic  response  to  atmospher¬ 
ic  forcing  is  constructed  to  assist  in  the  interpretation  of  the  obser¬ 
vational  results.  The  model  explores  the  consistency  of  seasona?  wind- 
stress  forcing  as  a  mechanism  for  creating  the  observed  seasonal  signals. 
The  model  reproduces  the  observations,  and  provides  a  means  for  extrac¬ 
ting  additional  information  from  the  observations.  The  windstress  model 
was  chosen  in  part  because  observations  of  the  large  scale  winds  are 
available.  An  alternate  model  might  hypothesize  that  the  observed  sea¬ 
sonal  signals  result  from  wintertime  convection  in  either  the  Labrador  Sea 
or  the  Mediterranean  Sea,  both  of  which  are  sources  of  intermediate  water 
in  the  Bay  of  Biscay.  It  is  equally  difficult  to  prove  or  disprove  that 
hypothesis,  since  the  amplitude  of  the  convection,  its  horizontal  struc¬ 
ture  and  extent  are  unknown;  neither  is  there  information  available  about 
the  expected  phase  relationship  between  the  convection  and  its  eventual 
influence  on  intermediate  waters  far  from  the  convective  source.  The 
short  horizontal  scale  of  the  salinity  anomaly  signal  implied  by  the 
observations,  and  confirmed  by  further  reduction  of  the  data  in  Chapter  3, 
argues  against  convection  as  a  source,  particularly  if  diffusion  is  im¬ 
portant.  However,  lacking  further  observations,  the  hypothesis  of  a 
convective  source  remains  untested. 
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Chapter  3 

3.1  Introduction 

The  observation  of  seasonal  signals  in  density  and  in  salinity 
anomaly  below  the  main  thermocline  introduces  several  questions  which 
cannot  be  answered  by  the  observations  alone :  What  are  the  horizontal 
scales  of  these  signals?  What  forces  them?  Why  do  they  appear  to  be 
stronger  below  the  thermocline?  Should  corresponding  signals  in  vertical 
displacement  of  density  surfaces  and  potential  energy  be  expected?  In 
this  chapter  an  attempt  to  answer  these  questions  is  made  through  the 
construction  of  an  analytical  model  of  atmospheric  (winds tress)  forcing 
and  oceanic  response.  The  model  is  highly  simplified,  in  order  to  be 
analytically  tractable,  but  by  hypothesis  contains  the  essential  physics 
of  the  problem.  (That  this  is  in  fact  true  must  be  verified  by  comparison 
of  the  neglected  terms  with  those  retained.)  The  model,  which  follows 
that  of  Bryan  and  Ripa,  1978,  assumes  that  the  forcing  consists  of  a  wave¬ 
form  which  may  be  propagating  eastward  or  westward,  or  standing,  and  whose 
wavelength,  phase,  amplitude  and  frequency  are  parameters.  The  forcing 
acts  upon  the  ocean,  through  a  mixed  layer  of  finite  depth,  as  a  body 
force  rather  than  a  stress.  The  model  ocean  is  continuously  stratified 
with  depth  variable  buoyancy  frequency  N  below  a  surface  mixed  layer; 
the  Coriolis  parameter  is  allowed  to  vary  with  latitude  in  the  usual 
3-plane  approximation.  There  _s  a  meridional  barrier  at  the  eastern 
edge  of  the  basin,  and  no  bottom  topography,  friction  or  diffusion.  The 
stratification,  mixed  layer  depth  and  total  depth  are  parameters.  The 
oceanic  response  to  the  forcing  consists  of  two  parts:  a  forced  wave  of 
wavelength  and  frequency  corresponding  to  the  atmospheric  forcing,  and 
a  sum  of  free  waves  required  to  satisfy  the  constraint  of  no  flow  into 
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the  boundary.  The  free  waves  are  Rossby  waves  with  meridional  wavenumber 
equal  to  that  of  the  forced  wave,  and  zonal  wavenumber  determined  by  the 
Rossby  wave  dispersion  relation.  The  forcing  scale  at  seasonal  period 
is  not  well  known,  but  it  is  certainly  as  large  as  the  basin.  Observat¬ 
ional  evidence  of  annual  windstress  forcing  is  discussed  in  section  4 
of  this  chapter.  The  free  response  zonal  scales  are  determined  primarily 
by  the  stratification,  which  is  well  known.  Once  the  dominant  free  wave 
scales  have  been  established,  the  model  can  then  predict  the  phase,  hori¬ 
zontal  and  vertical  displacements,  the  kinetic  and  potential  energies  of 
the  system. 

The  predicted  horizontal  displacements,  acting  on  a  tracer  field 
with  mean  horizontal  gradients,  will  deform  the  initial  tracer  field 
(reversibly,  in  the  absence  of  friction  and  diffusion)  so  that  the  progress 
of  the  disturbance  through  the  fluid  is  traced  by  the  anomaly  in  hori¬ 
zontal  structure.  Using  the  observed  mean  gradients  of  salinity  anomaly 
(As) ,  the  model  predicts  the  amplitude  and  phase  of  the  signal  in  As.  In 
order  to  compare  the  result  with  the  observed  salinity  anomaly  in  section  4, 


the  data  are  fit  to  a  progressive  wave  of  set  wavelength  and  frequency, 
and  the  amplitude  and  phase  calculated.  In  terms  of  particle  velocities, 
the  dominant  oceanic  response  away  from  the  coast  is  the  first  free  baro- 
clinic  wave,  so  only  that  wave  is  used  in  the  fit  to  the  observations.  A 
similar  analysis  is  performed  using  the  observed  vertical  displacements. 


Iiujt  H  t|  h  i.fc'l  n>  h  Ji  ill  i 
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3 . 2  Model 

The  governing  momentum  equations  are: 

(3.1a) 

(3.1b) 

with  x  and  y  positive  eastward  and  northward  co-ordinates,  u  and  v 
the  eastward  and  northward  velocities,  p  the  pressure  and  X  and  Y 
body  forces  which  parameterize  the  effect  of  the  windstress  through  the 
mixed  layer.  The  fluid  is  assumed  to  be  hydrostatic;  so  that 

1^-  =  -  pg  (3.1c) 


9u 

3t 


-fv  =  - 


—  ^  +  X 

p  3x  X 


3v 

3t 


+  fu 
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with  z  the  vertical  co-ordinate,  which  decreases  with  depth  from  a 
value  of  zero  at  the  surface.  It  is  also  assumed  that  the  fluid  is  in¬ 
compressible,  and  that  density  is  conserved  following  fluid  parcels.  The 
equation  of  continuity  is  then 


3u  .  3v  . 
3x  3y 


(3. Id) 


The  density  p  is  presumed  to  be  constant  through  the  mixed  layer  but 
a  function  of  x,  y,  z  and  t  below: 


p  =  p(z)  +  p^ (x,y , t) 


p  =  p  =  constant 
o 


0  >  z  > 


-h 
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P  = 


O  2 

p  +  —  N 

O  yg  O 


le 


Yh 


-eYZ] 
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where  h  is  the  depth  of  the  mixed  layer  and  H  is  the  depth  of  the 
fluid.  This  distribution  of  density  results  in  a  simple  buoyancy  fre¬ 
quency  profiles 
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p  dz 
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which  is  analytically  tractable.  The  density  equation  may  then  be  written 
as : 


Bp  ?QN2 
Bt  g 


W  =  0 


(3.1e) 


where  the  horizontal  advection  of  density  has  been  neglected;  w  is  the 
vertical  velocity. 

A  vorticity  equation  may  be  derived  by  subtracting  the  derivative 


with  respect  to  y  of  (3.1a)  from  the  derivative  with  respect  to  x  of 


(3.1b). 
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From  the  density  equation  (3.4),  and  using  (3.1c),  the  following  expression 
for  the  vertical  velocity  may  be  derived: 


3w  _  _B_  r_3_  1  &£ 
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(3.3) 
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If  the  geostrophic  streamfunction  p  is  defined  by 


p  fv  = 
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1 £ 
ax 
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then  (3.2)  may  be  written  in  terms  of  p: 


—  [V2  - 
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(3.4) 


with  V2  the  horizontal  laplacian  operator 
H 

constant. 

Equation  (3.4)  holds  provided 
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-A2  p  =  0, 


(3.5) 


which  allows  the  separation  of  the  solution  p  into  vertical  and  hori¬ 
zontal  structure  functions.  Thus  solutions  of  (3.4)  and  (3.5)  are  sought 
in  the  form: 


P  =  P(z)  Q(x)y) . 

(3.6a) 
(3.6b) 
(3.6c) 


The  boundary  conditions  imposed  are 
u  =  0  at  x  =  0 
w  =  0  at  z  =  0,  -H 
w,  p  continuous  across  z  =  -h. 


The  windstress  forcing  is  assumed  to  have  the  form 


X,  Y  =  (X  sin  Zy,  Y  cos  Zy)  exp  (:.(k£x  -  tot)). 


The  response  may  be  expressed  as  a  sum  of  a  forced  wave  and  an  infinite 
number  of  free  waves: 


p  =  Of  pf(z)  cos  Zy  exp  [i(k^x  -  tot))  + 

.  (3.7) 

Za  f>  ( z )  cos  Zy  exp[i(k  -  tot)] 

.  n  n  n 

n=0 

The  amplitudes  are  determined  by  application  of  the  boundary  con¬ 

dition  at  x  =  0.  The  vertical  structure  function  pf(z)  of  the  forced 
wave  is  determined  from  (3.5)  with  X  given  by  the  forcing  parameters 
according  to  (3.4): 


2  Bk1 
~  -k£  + 


The  structure  functions  6  (z)  of  the  free  waves  are  determined  from  (3.5) 

n 

for  discrete  values  of  X2  which  satisfy  the  constraint  of  continuous 

n 

pressure  across  the  base  of  the  mixed  layer. 

The  vertical  equation  (3.5)  can  be  rewritten  in  terms  of  the  trans¬ 
formed  z-variable  s  =  exp(yz/2): 


2  a  a>  2*,2_  n 

sp  -sp  +  s.pp  =  0 


-h  >  z  >  -H 
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which  can  be  solved  in  terms  of  first  order  Bessel  functions.  Since 
the  boundary  conditions  on  the  vertical  equation  are  easily  expressed 
in  terms  of  vertical  velocity,  it  is  simplest  to  relate  pressure  to 
vertical  velocity  w  and  solve  the  problem  in  w.  In  doing  that  the 
barotropic  mode  is  excluded;  in  order  to  retain  both  the  barotropic  mode 
and  the  simple  boundary  conditions,  solutions  in  p  a-<d  in  w  will  be 
used,  and  a  consistent  conversion  between  them  established.  This  pro¬ 
cedure  departs  from  that  used  by  Bryan  and  Ripa,  who  ignore  the  baro¬ 
tropic  wave  altogether.  Pressure  and  vertical  velocity  are  related  by: 


or 


9w 9 3_  .  1  3p, 

3z  3t  3z  N2p  3z‘ 
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The  equation  for  w  is  then 
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with  solutions  of  the  form 
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j  and  Y  Bessel  functions, 
o  o 

In  addition  to  boundary  conditions  at  the  top  and  bottom,  the  ver¬ 
tical  velocity  and  its  first  derivative  are  constrained  to  be  continuous 
across  z  =  -h.  The  mixed  layer  solution  for  ^  is  simply  proportional 
to  z.  Applying  the  constraint  on  w  at  z  =  -h,  the  solution  for 


0  >  z  >  -h 


w  =  -  z  c  t?  (y  s*)/h 
n  n  ro  rn 

w  =  c  t>  (y  s)  h  >  z  >  -H  (3.10) 

n  n  r  o  n 

s*  =  s(-h)  =  exp[-yh/2] 


The  condition  w  =  0  at  z  =  -H  is  satisfied  separately  by  each  mode  if 
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e  =  s(-H)  =  exp  [-yH/23 


(3.11) 


The  eigenvalues  y^  are  determined  by  application  of  the  constraint  on 
r\ 

w  at  z  =  -h: 
z 


y  s* 
Hn 


(3.12) 


The  condition  of  continuous  across  z  =  -h  is  equivalent  to  con¬ 

tinuous  p.  Using  (3.8)  the  solution  in  terms  of  p  is: 


0  >  z  >  -h 
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The  eigenfunctions  p^  are  orthonormal,  and  the  normalizing  constant 

determines  c  ;  that  is, 
n 
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(with  6  the  Kroenecker  delta)  implies  that 
mn 
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(3.14a) 


for  n  >  1,  and 
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(3.14b) 


For  the  forced  wave,  if  the  wave  is  eastward-moving  A*  and  there¬ 
fore  y*  are  negative  and  the  solution  for  w^  is  in  terms  of  modified 
Bessel  functions 


wf  =  -zcf  (Hfs*)/h 
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For  a  westward  moving  forced  wave  the  solution  is 


0  >  z  >  -h 
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(3.15) 


wf  =  -z  cf  If  Q  (pfs*)/h 


0  >  z  >  -h 


The  constants  bf  and  b_  are  determined  by  requiring  that  w 
vanish  at  z  =  -H.  The  coefficients  c^  and  c^  are  determined  by 
assuming  that  the  vertical  velocity  at  the  base  of  the  mixed  layer  due 
to  the  wind  is  given  by  the  divergence  of  the  Ekman  mass  flux  within 
the  mixed  layer: 
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The  coefficients  cf  are  further  adjusted  so  that  the  amplitudes  of 
wf  and  wf  are  equal  at  the  base  of  the  mixed  layer.  Thus 
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The  forced  wave  streamfunctions  corresponding  to  (3.17)  are: 


A 


'V1 


(3 


where 


+  ei”  V- 


106 


Finally,  in  order  to  determine  the  amplitudes  a  of  the  free  waves 

n 

the  boundary  condition  at  x  =  0  is  applied: 


u(x  =  0)  =  0 


or 


Z  sin£y  p^(z) 


M 

Zl  sinJly  a  j?  (z). 
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Applying  the  orthogonality  condition  (3.13): 


e 


which  can  be  evaluated  analytically  (Abramowitz  and  Stegun,  1972,  p.  484). 
The  solution  depends  upon  the  direction  of  the  forcing;  for  the  eastward 
propagating  wave: 
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(where  subscript  zero  refers  to  the  barotropic  mode)  while  for  the  west¬ 
ward  propagating  wave 
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The  boundary  condition  of  no  zonal  velocity  at  x  =  0  assures  that  there 
is  no  energy  flux  into  the  boundary.  As  x  approaches  -<*>,  the  energy 
from  the  free  waves,  which  exist  only  because  of  the  presence  of  the  meri¬ 
dional  barrier,  roust  be  propagating  away  from  the  barrier,  or  westward. 
This  radiation  condition  requires  that  only  the  free  waves  with  westward 
group  velocity  be  included  in  the  solution.  The  dispersion  relation  is 

a  quadratic  in  k,  so  that  two  solutions  are  possible  for  each  X  .  The 

n 

shorter  wave  has  eastward  group  velocity,  the  longer  wave  westward  group 
velocity.  The  long  waves  are  therefore  chosen  for  the  solution. 
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3.3  Model  results 

The  parameters  used  in  the  calculation  are  listed  in  Table  3.1 
Discussion  of  the  windstress  observations  which  were  used  in  choosing 
the  forcing  parameters  is  deferred  to  section  4.  The  mixed  layer  depth, 
total  depth  and  the  stratification  parameters  are  taken  from  the  Bay  of 
Biscay  observations.  The  eigenvalues  of  the  first  four  modes  were  cal¬ 
culated  using  a  Newton- Raph son  iteration  of  (3.12);  they  are  tabulated 
in  Table  3.2.  For  .n=2  and  higher,  the  modes  are  trapped  within  20  km 
of  the  coast.  (Recall  that  the  observations  extend  roughly  700  km  from 
the  coast.)  The  dispersion  relation  is: 


k 

n 


-  4  (A2  +  fc2) 
n 


For  modes  with  n  >  2  has  two  complex  roots;  the  root  which  decays 
away  from  the  coast  must  be  chosen.  Because  the  higher  modes  are  strongly 
trapped  only  the  barotropic  and  first  baroclinic  modes  are  included  in 
the  calculations  which  follow.  The  wavelengths  and  phase  speeds  of  these 
modes  are  also  given  in  Table  3.2.  The  first  baroclinic  wave  has  a  wave¬ 
length  of  391  km.  The  barotropic  wave  is  much  longer  than  the  basin  and 
therefore  has  negligible  displacements.  The  vertical  and  horizontal 
displacements  associated  with  the  first  baroclinic  wave  are  much  larger 
than  those  associated  with  the  forced  wave  and  the  barotropic  wave. 

As  a  result  the  effective  difference  in  terms  of  displacements  between 
propagating  and  standing  windstress  forcing  is  the  amplitude  of  the 
free  response:  for  a  standing  wave  the  location  of  the  nodes  of 
the  wave  pattern  are  important  in  determining  the  response  amplitude. 

Since  the  horizontal  structure  of  the  annual  windstress  is  not  well 
known,  the  eastward  propagating  wave  was  chosen  as  representative  and  the 
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TABLE  3.1 
Model  Parameters 


Parameter 

Description 

Value 

h 

Mixed  layer  depth 

400  ra 

H 

* 

Total  depth 

5000  m 

2ir/N0 

Buoyancy  period  at  h 

1.35  x  ios  sec 

1/Y 

• 

e-folding  depth  for  stratification 

1300  m 

2TT/kf 

• 

Zonal  wavelength  of  forcing 

12000  km 

2ir/l 

: 

Meridional  wavelength  of  forcing 

6000  km 

2tt/w 

: 

period 

1  year 

B 

Beta 

1.6  x  io  11  (m  sec' 

f 

: 

Coriolis  parameter 

10  4  sec  1 

• 

Meridional  wind  stress 

10.  x  10  2  pascal 

TX 

« 

• 

Zonal  windstress 

—  4.x  10  2  pascal 
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TABLE  3.2 


Eigenvalues  and  Corresponding  Wavelengths  and  Phase  Speeds 


2n/k 

(x)/k 

MODE 

^n 

km 

m  sec  1 

forced 

.647 

12000. 

.381 

Barotropic  Mode 

.049 

-4.1  x  105 

-13.0 

Baroclinic  Modes 

1 

3.59 

-391. 

-.0125 

2 

7.57 

3 

11.7 

4 


16.0 


Ill 


results  which  follow  are  for  that  case. 

The  streamfunctions  p^,  pq,  3^  are  plotted  in  Fig.  3.1,  illus¬ 
trating  the  relative  amplitudes  of  the  barotropic  and  baroclinic  components. 
Recall  that  the  boundary  condition  at  x  =  0  requires  that  the  total 
streamfunction  vanish  there.  As  a  result,  no  signal  in  sea  level  at  the 
coast  is  expected. 

Displacement  amplitudes  £  as  a  function  of  depth  were  calculated 
as 

✓s 

§§■  =  fc(z) 

at. 

or  (3.21) 


and  are  plotted  in  Fig.  3.1  for  both  the  forced  and  the  first  free  inodes. 
The  meridional  velocity  amplitude 

v  =  k  p(z)  (3.22) 

is  plotted  in  Fig.  3.3,  for  the  forced  and  first  baroclinic  modes.  The 
meridional  velocity;  of  the  barotropic  mode  is  two  orders  of  magnitude 
smaller  them  that  of  the  forced  mode. 

Since  the  free  wave  zonal  velocity  is  much  smaller  than  the  meri¬ 
dional  velocity,  the  salt  conservation  equation  with  no  diffusion  is 
approximately 


or 


9s 

at 


+  v(x) 


(3.23) 
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S(x,y,z,t.)  =  So(x,n,0)  •  S  (z) 


n  -  y 


(3.24) 


For  a  simple  initial  field  Sq,  With  salinity  anomaly  a  linear  function 
of  y  and  an  arbitrary  function  of  x: 


S(x,y,z,t)  =  (So(x,y,0)  +  Sy  *  j  vdt]  S(z). 


Using  the  observed  values  of  from  the  Bay  of  Biscay  dara, 

one  may  obtain  an  estimate  of  the  salinity  anomaly  variation  amplitude: 


&  .  -<s  )  „  % 

y  obs  no 


(3.25) 


which  is  plotted  in  Fig.  3.6  for  those  levels  where  (Sy^0ks  -*-s 
available . 

Note  that  certain  phase  relationships  predicted  by  the  model  may 
be  tested.  If  the  windstress  varies  as 


t  =  cos  $  +  i  sin  $ 


$  =  kfx  -  tot  -  <j> 


then  the  displacement  and  the  salinity  anomaly  fields  will  vary  as 
(cos  $  +  i  sin  <J>) ,  whereas  the  velocity  will  be  90°  out  of  phase  with  the 
forcing.  The  observed  phase  4>  in  salinity  anomaly  and  displacement  can 
be  compared  to  the  observed  phase  of  the  annual  windstress#  as  discussed 
in  the  following  section. 
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3.4  Comparison  of  Model  and  Observations 

There  are  two  concerns  in  this  section:  a  discussion  of  the  ob¬ 
servation  of  windstress  in  the  North  Atlantic  which  were  used  to  de¬ 
termine  the  choice  of  forcing  parameters  in  the  model  and  the  phase  of 
the  forcing;  and  a  further  reduction  of  the  oceanic  observations  of 
Chapter  2  which  allows  direct  comparison  of  the  model  predictions  and 
the  observational  results. 

There  exists  a  comprehensive  meteorological  data  set  for  the 

North  Atlantic  from  1948  to  1972  (Bunker  and  Worthington,  1976) .  Charts 

of  the  average  annual  windstress  over  the  North  Atlantic  were  prepared 

by  Bunker  and  Worthington;  to  estimate  the  amplitudes  and  wavelengths 

of  the  annual  windstress  forcing  appropriate  to  this  problem  zonal  and 

x  y 

meridioru  >.  sections  of  x  and  x  (along  45°N  and  10°W)  were  taken 
from  those  charts  (see  Fig.  3.4).  The  entire  25  year  record  of  windstress 
was  examined  for  evidence  of  pjiase  propagation,  with  no  conclusive  results. 
There  is  some  evidence  that  the  zonal  windstress  maximum  shifts  eastward 
in  summer  and  fall  relative  to  its  position  in  winter  and  spring  (Fig.  3.5). 
The  forcing  wavelengths  used  in  the  model  were  estimated  from  Fig.  3.5. 

The  model  results  are  not  sensitive  to  the  choice  of  forcing  scales  within 
reasonable  limits.  The  forcing  amplitudes  were  chosen  such  that  the 
model  As  response  amplitude  was  approximately  that  given  by  the  obsei- 
vations,  with  the  ratio  of  t  /  t  as  given  by  the  wind  observations  held 
constant.  The  long  term  average  amplitude  of  seasonal  change  in  windstress 
estimated  from  the  difference  between  winter  and  summer  25  year  means  is 

roughly  5.  x  10  2  pascals  (.5  dyne  cm  2)  for  tX  and  2  *  10  2  pascals 

v  x 

for  x*  along  45°N  between  10°  and  70°W.  One  standard  deviation  of  x 


120 


Figure  3.4  (a)  Meridional  sections  along  10°W  and 

(b)  Zonal  sections  along  45°N  of  the  annual  average 

x  y 

windstress  components  x  and  x  .  Values  taken  from 
Bunker  and  Worthington  (1976) . 


Figure  3.5  Seasonal  averages  using  all  the  Bunker  and  Worthington  (1976) 


data  (1948-1972)  of  TX  along  45°N,  by  10°  Marsden  square. 


— 2  y  “2 

is  13.  x  io  pascals  and  of  t  9.7  x  10  pascals  for  the  winter 
season  in  the  ten-degree  Marsden  Square  corresponding  to  40°  to  50 °N, 

0°  to  10°W.  The  standard  deviations  there  for  the  summer  season  are 
5.2  x  io  2  pascals  and  3.1  x  10  2  pascals  for  tX  and  respec¬ 

tively.  Similar  variances  are  found  along  45°N  between  10°  and  70°W. 

The  windstress  amplitudes  used  in  the  model  are  10  x  10  2  and  -4  x  10  2 
x  y 

pascals  for  t  and  t  ,  within  a  standard  deviation  of  the  25-year 
means. 

The  phase  of  the  seasonal  cycle  of  windstress  was  calculated  by 
Firing  (1978)  for  different  longitude  bands.  At  38  °N,  between  30°W  and 
40°W,  he  calculated  that  the  maximum  in  windstress  occurred  55  days  from 
January. 

The  model  predicts  westward  propagating  free  disturbances  of  domi¬ 
nant  wavelength  391  km,  frequency  1  cpy  and  phase  corresponding  to  a 
maximum  at  55  days  from  1  January.  In  Chapter  2  the  data  were  fitted  to 
polynomials  in  x  and  y  to  study  the  horizontal  structure  of  salinity 
anomaly  and  vertical  displacement.  That  procedure  is  inappropriate  for 
comparing  the  model  and  the  observations,  since  the  propagation  of  waves 
through  the  field  results  in  a  complex  interaction  of  space  and  time  infor¬ 
mation  which  cannot  be  sorted  out  by  fitting  to  polynomials.  Instead, 
the  data  were  fitted  to  traveling  plane  waves  of  the  form 

cos  (kx  -  wt  -  <*>)  • 

In  that  analysis,  the  steady  horizontal  structure  of  the  salinity 
anomaly  field  was  removed.  The  residuals  were  then  fitted  to 

5  =  3^  sin^x  -  ut)  +  a2  cos (k^x  -  cat)  (3.26) 


with 
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k1  =  27T/391  km 

a)  =  2tt/365  days 

a,  =  -a  sin  d> 

Is 

a„  =  a  cos  d> 

2  s 

and  x  measured  as  actual  horizontal  distance  from  the  600  m  isobath 
along  the  eastern  boundary.  The  600  m  isobath  does  not  correspond  to  a 
meridional  barrier;  the  implications  of  that  choice  for  x  =  0  are  dis¬ 
cussed  in  section  3.5. 

The  regression  was  done  using  all  cruises  at  each  level/  and  sep¬ 
arately  by  cruise.  In  order  to  determine  the  amplitudes  with  statistical 
confidence  greater  than  90  per  cent,  the  best  estimate  of  phase  was  used 
for  <J>,  and  the  amplitude  alone  was  fitted,  with  fixed  phase. 

The  best  estimate  of  <j>  corresponds  to  a  maximum  at  71  days  from 
1  January.  The  estimated  erro$  in  the  phase  determination  is  ±.3  rad 
(20  days) ,  so  that  the  calculated  phase  agrees  with  that  of  Firing  within 
the  errors  of  phase  determination.  The  overall  amplitude  a^  as  a  func¬ 
tion  of  z  is  plotted  in  Fig.  3.6,  along  with  the  amplitudes  for  cruises 
which  show  strong  signals,  and  the  model  prediction  for  an  eastward 

forced  wave.  Most  of  the  variance  in  a  is  found  in  four  cruises: 

s 

Phygas  23,  31,  42,  44;  the  others  contribute  very  little  signal.  The 

lack  of  signal  may  be  due  in  part  to  insufficient  sampling,  as  there  are 

significantly  fewer  stations  in  those  cruises  (see  Table  2.1).  For 

Phygas  42  the  amplitude  a  is  positive,  although  it  is  plotted  as  -a 

s  s 

in  Fig.  3.7.  The  vertical  structure  is  similar  in  the  four  cruises  with 
large  signals,  and  compares  very  well  with  the  model  prediction. 


;  »  yjqm  ggg^gTr;^ 


Figure  3.6  Salinity  anomaly  amplitudes:  heavy  solid  line  is  the  model 


prediction,  dased  line  is  the  overall  value  (excluding 


cruises  42  and  43)  and  the  remaining  thin  lines  are  cruises 


24,31,44,  and  the  negative  of  42,  as  noted. 
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In  order  to  determine  how  sensitive  the  regression  fit  is  to  the 
chosen  parameters  k,  to,  and  <{),  tests  were  run  using  a  range  of 
values  for  each  parameter,  holding  the  other  two  at  their  chosen  values. 
The  ranges  covered  were 

2iv/k  (391  km)  :  300  to  900  km,  steps  of  100  and  200  km 

2ir/<p  (12  months):  10  to  14  months,  steps  of  1  month 
(J)  -  365/2  it  (71  days):  50  to  130  days,  steps  of  20  days 
In  all  cases  the  ratio  of  the  amplitude  to  its  regression  standard  dev¬ 
iation  was  a  maximum  for  the  chosen  parameters,  and,  with  the  exception 
of  <{>  corresponding  to  50  to  90  days,  none  other  than  the  one  with  the 
three  chosen  parameters  were  significant  at  the  95  per  cent  confidence 
level  (see  Fig.  3.7).  The  estimated  error  in  the  determination  of  <{) 
given  above  is  taken  from  these  results. 

A  similar  analysis  was  carried  out  for  the  displacement ■ field, 
although  there  were  no  identifiable  mean  horizontal  gradients  to  be  re¬ 
moved.  The  observed  displacement  amplitudes  (z)  are  plotted  in 
Fig.  3.8  along  with  the  model  prediction  and  values  for  Phygas  31  md  24 
which  showed  the  most  consistent  single  cruise  results. 

The  model  predicts  little  seasonal  change  in  APE  per  unit  mass  which 
is  computed  as  an  average  over  the  area.  Since  the  zonal  extent  of  the 
d.ca  is  longer  than  the  zonal  wavelength  of  dominant  interior  response, 
the  predicted  APE  varies  little  with  time.  The  model  APE  per  unit  mass 
does  have  vertical  structure  similar  to .that  observed,  with  a  maximum 
value  at  about  1400  db.  The  amplitude  of  the  model  APE  is  more  than  an 
order  of  magnitude  smaller  than  observed;  this  is  not  surprising,  since 
the  root  mean  square  displacements,  calculated  in  Chapter  2,  which  include 
non-seasonal  as  well  as  seasonal  variability,  are  5  to  10  times  as  large 
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Figure  3.8  Displacement  amplitudes:  heavy  solid  line  is  the  model 
prediction,  dashed  line  is  the  overall  value  (excluding 
cruises  42  and  43)  and  the  thin  lines  are  cruises  31  and  24. 


DISPLACEMENT  AMPLITUDES,  m 


as  the  annual  displacements  a^  calculated  above.  (See  Figures  3.9 
and  2.5.)  An  estimate  of  the  observed  annual  component  of  APE  can  be 
made  using  the  displacement  amplitudes  a^: 


APE 


annual 


=  .5  N2  'aT 
TT 


(3.27) 


which  is  in  agreement  with  the  model  prediction  (Fig.  3.9) 

Of  the  motions  predicted  by  the  model,  only  the  forced  wave  has 
sufficiently  long  horizontal  scale  to  explain  the  changes  in  the  reference 
steric  field,  which  is  independent  of  horizontal  position.  The  observed 
seasonal  signal  in  the  reference  steric  field  corresponding  to  a  vertical 
displacement  of  0(-15  m)  at  the  thermocline  is  not  well  explained  by  the 


model,  which  predicts  a  displacement  of  0(+5  m)  at  the  base  of  the  mixed 
layer,  decreasing  monotonically  below,  for  the  forced  wave. 


Figure  3.9  APE  per  unit  mass:  heavy  solid  line  is  the  model 


prediction  (.5  N2  £2)  and  triangles  are  APEannua^  given 


by  equation  (3.27). 


id  - 


?; 


136 


3.5  Discussion 

In  this  chapter  an  analytical  model  of  seasonal  atmospheric  forcing 
and  oceanic  response  was  constructed  in  an  attempt  to  answer  a  number  of 
questions  raised  by  the  observational  results  of  Chapter  2.  Foremost 
among  those  questions  was  that  of  the  horizontal  scale  of  the  seasonal 
signal  in  salinity  anomaly.  Secondly,  the  question  of  forcing  was  posed: 
what  types  of  atmospheric  forcing  might  result  in  the  observed  seasonal 
signals,  with  regard  to  both  horizontal  scale  and  amplitude?  Finally, 
what  other  variables  might  be  expected  to  show  corresponding  signals  con¬ 
sistent  with  the  model  forcing? 

In  the  model  ocean,  the  response  consists  of  two  parts:  a  forced 
wave  with  scales  comparable  to  the  atmospheric  forcing,  and  a  set  of 
free  waves  required  uo  satisfy  the  boundary  condition  at  a  meridional 
eastern  barrier.  The  dominant  free  response  away  from  the  coast  was 
found  to  be  the  first  baroclinic  mode  of  wavelength  390  km.  The  baro- 
troplc  wave  is  longer  than  the  basin  width,  and  has  very  small  displace¬ 
ments.  The  second  and  higher  modes  are  found  to  be  trapped  within  20  km 
of  the  coast.  Thus,  without  any  assumptions  about  the  scale  or  amplitude 
of  the  forcing,  the  scale  of  the  free  annual  signal  in  the  interior  has 
been  established,  and  it  has  been  determined  that  no  shorter  Rossby  waves 
can  escape  the  coastal  boundary  region. 

In  order  to  determine  the  relative  importance  of  the  forced  and 
free  waves  in  the  interior  further  assumptions  about  the  forcing  r<ust 
be  made.  For  simplicity,  the  model  forcing  is  assumed  to  be  windstress 
forcing;  this  is  mainly  because  more  data  on  the  distribution  of  wind- 
stress  are  available  than  on  other  types  of  atmospheric  forcing.  However, 
since  the  windstress  forcing  is  further  assumed  to  act  as  a  body  force 


137 


through  the  mixed  layer,  so  that  it  affects  the  remainder  of  the  water 
column  only  through  an  imposed  vertical  velocity  at  the  base  of  the  mixed 
layer,  the  model  is  certainly  compatible  with  other  types  of  forcing. 

The  horizontal  scales  of  the  model  forcing,  estimated  from  the  observations 
compiled  by  Bunker  and  Worthington  (1976),  correspond  to  wavelengths  of 
12,000  and  6,000  km  in  the  zonal  and  meridional  directions  respectively. 

The  amplitudes  used  in  the  model  are  .1  and  -.04  pascals  (1.  and  -0.4 
dyne-cm  )  for  t  and  T  respectively.  The  application  of  this 
forcing  to  the  model  ocean  results  in  estimates  of  vertical  and  horizontal 
displacements  for  the  free  and  forced  waves.  The  meridional  particle 
speed  of  the  first  baroclinic  wave  is  approximately  equal  to  the  phase 
speed  (0(.01  m  sec  *);  for  the  barotropic  wave  the  meridional  velocity  is 
OdO-**  m  sec  x)  and  for  the  forced  wave  0(10  3  m  sec  *). 

The  change  in  the  zonal  structure  of  the  salinity  anomaly  field 
predicted  by  the  model  is  estimated  by  the  product  of  the  model  horizontal 
displacement  and  the  observed  steady  meridional  gradient  in  salinity 
anomaly.  For  comparison  with  the  model  prediction,  the  observations  were 
fit  to  a  traveling  plane  wave  of  wavelength  390  km.  The  phase  and  ampli¬ 
tude  were  both  estimated  from  the  observations.  The  estimated  phase  in 
salinity  anomaly  agrees  with  the  observed  phase  of  windstress  ..ithin  the 
errors  of  phase  determination.  The  vertical  structure  of  the  model  and 
the  observed  salinity  anomaly  amplitude  agree  very  well.  The  amplitude 
of  the  signal  using  all  cruises  except  Phygas  42  is  a  maximum  of  -12  pprr 
compared  to  -16  ppm  predicted  by  the  model.  The  individual  cruises  24, 
31,  42,  44  have  larger  absolute  amplitudes,  of  order  25  to  30  ppm.  Phygas 
42  is  the  only  cruise  for  which  the  amplitude  is  positive;  its  vertical 
structure  is  indistinguishable  from  the  other  cruises.  A  positive 
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amplitude  is  predicted  by  the  model  only  for  positive  windstress  curl. 
Although  not  generally  the  case  in  this  part  of  the  North  Atlantic,  posi¬ 
tive  curl  represents  something  like  one  to  two  standard  deviations  of  the 
observed  stresses  from  the  25  year  means,  at  the  model  horizontal  scales. 

The  observed  vertical  displacements  were  also  fit  to  a  traveling 
plane  wave  of  390  km  wavelength;  the  amplitude  and  vertical  structure  of 
the  observed  signal  agrees  well  with  that  of  the  model.  This  is  an  in¬ 
dependent  test  of  the  consistency  of  the  model  and  of  the  forcing  ampli¬ 
tude.  The  observed  vertical  displacement  amplitude  is  0(15  m)  for  the 
390  km  wave. 

The  observations  of  seasonal  variability  in  the  zonal  structure  of 
salinity  anomaly  and  vertical  displacement  are  therefore  consistent  with 
the  proposed  model  of  windstress  forcing,  with  respect  to  amplitude, 
phase,  vertical  and  horizontal  structure.  The  observed  seasonal  signal 
in  the  reference  steric  field  and  the  apparently  seasonal  oscillations 
in  the  APE  below  the  main  thermocline  are  not  adequately  explained  by  the 
model;  in  the  former  case  the  observed  signal  is  a  factor  of  3  too  small 
and  has  the  wrong  phase,  and  in  the  latter  the  amplitude  of  the  observed 
signal  is  more  than  an  order  of  magnitude  larger  than  the  model 
prediction. 

The  model  presented  in  this  chapter  is  simplified  in  order  that  it 
be  analytically  tractable.  Specifically,  friction,  diffusion,  and  bottom 
topography  are  ignored,  the  complex  geometry  of  the  coastline  has  not  been 
incorporated,  and  the  region  is  assumed  to  be  semi-infinite,  with  no 
western  boundary.  The  neglect  of  bottom  topography  is  not  serious,  as 
most  of  the  observations  are  from  the  area  of  the  Biscay  Abyssal  Plain. 
Estimates  of  the  size  of  the  friction  term  relative  to  the  time  derivative 
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term  in  the  momentum  equation  are  found  in  Table  3.3.  Friction  is  ap¬ 
parently  not  important  for  the  first  baroclinic  mode.  Lateral  diffusion 
has  a  stronger  effect  on  the  solution:  from  the  estimate  given  in 
Table  3.3  one  might  expect  an  error  of  25%  in  the  predicted  salinity 

anomaly  signal  due  to  the  neglect  of  diffusion,  for  a  diffusion  coeffi- 
3  2  -1 

cient  of  10  m  sec  .  For  the  present  case  that  error  is  a  maximum  of 
7  ppm^less  than  the  measurement  error  in  salinity. 

The  coastline  geometry  used  in  the  model  is  a  simple  meridional 
barrier  which  is  a  reasonable  approximation  of  the  eastern  coast  of  the 
North  Atlantic  over  a  scale  of  5000  km.  Locally,  however,  the  coastal 
geometry  is  more  complicated;  in  the  Bay  of  Biscay  it  will  be  noted  that 
the  eastern  boundary  (taken  as  the  shelf  break)  lies  at  an  angle  of  per¬ 
haps  45°  west  of  north.  If  the  model  assumed  an  angled  barrier  of  in¬ 
finite  extent,  then  there  is  a  critical  angle  for  which  all  free  baroclinic 
waves  are  coastally  trapped  (Philander,  1978).  For  the  parameters  used  in 
this  model  that  critical  angle  is  30°  west  of  north.  However,  the  con¬ 
straint  of  infinite  extent  is  unrealistic,  especially  since  the  north- 
south  extent  of  the  angled  coastline  is  considerably  less  than  the 
north-south  scale  of  the  forcing.  The  observations  suggest  that  the  pri¬ 
mary  effect  of  the  angled  boundary  is  an  alteration  of  the  phase  of  the 
free  response,  since  the  best  fit  of  the  data  to  cos(kx  -  tot  -  <j>)  used 
the  distance  from  the  600  m  isobath  as  the  x~coordinate.  It  should  be 
noted  that  in  the  absence  of  a  model  which  includes  the  coastal  geometry, 
this  suggestion  is  not  justified  theoretically.  A  proper  solution  to  the 
theoretical  problem  might  involve  a  separation  of  the  region  into  a  tri¬ 
angular  section  bounded  by  the  Iberian  Peninsula  to  the  south,  the  coast 
of  France  to  the  east  and  north  and  the  meridian  at  9°W;  and  the  open 
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ocean  west  of  the  triangular  section.  Separate  solutions  for  each  area, 
with  matching  conditions  along  9°W,  could  then  be  found.  The  solution  in 
the  triangular  section  will  most  likely  involve  Bessel  functions  in  the 
horizontal  structure,  which  introduces  a  practical  difficulty  in  the  com¬ 
parison  of  such  a  model  with  data. 

The  model  assumes  chat  the  region  is  semi-infinite,  implying  that  the 
western  boundary  is  sufficiently  far  way  that  it  has  no  effect  on  the 
solution  locally.  In  the  absence  of  friction  both  waves  will  eventually 
reach  a  western  boundary.  A  small  amount  of  friction  would  make  the 
western  boundary  effectively  much  farther  away  for  the  relatively  short 
baroclinic  wave.  Even  if  both  waves  generated  at  the  eastern  boundary 
actually  encounter  the  western  boundary,  and  reflect,  the  reflected 
waves,  with  eastward  group  velocity,  will  have  such  small  scales  that 
only  a  small  amount  of  friction  would  prevent  them  from  propagating  back 
across  the  basin  to  the  eastern  barrier.  Therefore  the  assumption  of  a 
semi-infinite  region  is  acceptable,  provided  a  small  amount  of  friction 
is  included  in  the  argument. 

In  summary,  the  observations  indicate  seasonal  signals  in  salinity 
anomaly  of  amplitude  ±25  ppm  and  vertical  displacement  of  amplitude±15  m, 
with  a  zonal  wavelength  of  390  km  and  a  phase  of  71  days  from  1  January, 
with  statistical  confidence  exceeding  90%.  The  model  presented  in  this 
chapter,  for  which  the  forcing  is  an  annual  signal  in  windstress  with 
horizontal  wavelengths  of  12000  and  6000  km  in  the  zonal  and  meridional 
directions,  amplitudes  of  .1  and  -.04  pascal  and  phase  corresponding  to 
a  maximum  at  55  days  from  1  January,  predicts  the  observed  salinity 
anomaly  and  vertical  displacement  signals  with  regard  to  wavelength, 
vertical  structure,  amplitude  and  phase  within  the  errors  of  the 


observations. 
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CHAPTER  4 

This  work  had  two  primary  and  three  corollary  objectives.  The  primary 
objectives  were:  1)  to  demonstrate  observational  evidence  for  seasonal 
variability  below  the  main  thermocline,  and  2)  to  determine  theorol ically 
the  types  of  atmospheric  forcing  which  are  consistent  with  the  observa¬ 
tions.  The  corollary  objectives  were:  1)  to  identify  precisely  from  ob¬ 
servations  those  material  surfaces  which  are  appropriate  to  the  study  of 
tracers  in  the  ocean,  2)  to  establish  a  general  technique  for  examining 
the  horizontal  distribution  of  tracers  along  those  surfaces,  and  3)  to  de¬ 
termine  from  observations  the  oceanic  available  potential  energy  (APE)  with 
particular  emphasis  on  identifying  the  errors  in  the  estimate  of  APE.  In 
this  chapter  the  success  of  those  endeavors  is  analyzed,  the  implications 
of  the  results  are  discussed,  and  suggestions  are  made  for  possible  further 
work . 

4 . 1  Results 

The  hyd.-egr aphic  observations  in  the  Bay  of  Biscay  have  been  examined 
by  calculating  the  adiabatically  leveled  reference  steric  field  appropriate 
to  each  cruise,  the  distribution  of  the  tracer  salinity  anomaly  along  those 
reference  steric  surfaces,  the  vertical  displacement  of  the  observed  steric 
surface  from  the  reference  surface  for  each  cruise  as  a  function  of  depth 
and  horizontal  position,  and  the  available  potential  energy  (APE) .  The 
three-year  duration  of  the  data,  together  with  the  sampling  interval  of 
three  months  and  the  large  spatial  extent  of  the  data  in  both  horizontal  di¬ 
mensions  permits  the  evaluation  of  seasonal  variability  in  salinity  anomaly, 
vertical  displacement,  and  the  reference  steric  field  itself  with  high 
statistical  confidence.  The  salinity  anomaly  field  experiences  seasonal 
fluctuations  which  appear  as  westward  propagating  plane  waves  with  zonal 
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wavelength  3-0  (±50)  km.  The  disturbance  is  180°  out  of  phase  with  the 
seasonal  component  of  the  zonal  wind  stress  which  has  a  maximum  value  at 
55  days  from  1  January.  The  computed  phase  in  salinity  anomaly  corres¬ 
ponds  to  a  minimum  at  71  (±30)  days  from  1  January.  The  amplitude  of  the 
seasonal  signal  in  salinity  anomaly  varies  with  depth,  with  a  maximum 
amplitude  of  30  (±10)  ppm  between  1000  and  1200  db  (Fig.  3.6) .  There  is 
also  a  seasonal  signal  in  the  vertical  displacement  field,  with  the  same 
horizontal  structure  observed  in  salinity  anomaly.  The  signal  in  dis¬ 
placement  also  varies  with  depth,  with  structure  corresponding  to  that 
expected  in  vertical  displacement  for  a  first  baroclinic  mode.  The  max¬ 
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imum  amplitude  is  30  (±7)  db. 


The  reference  steric  field  oscillates  seasonally  with  an  amplitude  of 

1  x  10  8 ( ± .  5  x  10  8)m\g  1  (1  ±  .5  x  10  5cm8gm  “S  .  Since  the  reference 

steric  field  is  independent  of  horizontal  position,  it  appears  that  this 

fluctuation  occurs  over  horizontal  scales  at  least  as  large  as  the  area 

a 

covered  by  the  observations.  The  phase  of  that  oscillation  corresponds 

to  a  maximum  at  130  (±30)  days  from  1  January. 

The  APE  per  unit  mass  in  this  region  has  two  maxima  in  the  vertical: 

one  corresponding  to  the  thermocline,  and  the  second  at  about  1400  db.  The 

thermocline  APE  is  comparable  in  magnitude  to  that  computed  in  MODE  (Bray 

-4  -1  2  -2 

and  Fofonoff,  1980)  with  a  typical  values  of  70  x  10  J  kg  (70  cm  sec  ) 
(Fig.  2.6).  The  second  maximum  is  comparable  in  magnitude  to  the  thermo¬ 
cline  maximum  for  most  cruises  and  exceeds  the  thermocline  values  for  some 
cruises  (Fig.  2.9).  There  is  visual  similarity  in  the  vertical  structure 
for  cruises  which  occur  in  the  same  seasons;  Phygas  31  and  41  which  are 
winter  cruises  and  Phygas  43,  a  summer  cruise  have  the  strongest  deep  max¬ 
ima  relative  to  the  reference  level  at  1100  db,  with  amplitudes  of 

-4  -1 

approximately  30  x  10  J  kg  relative  to  the  reference  level  APE. 
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Measurements  of  eddy  kinetic  energy  from  moored  current  meters  located 
at  47°N,  10°W  allows  comparison  of  the  ratio  of  APE  to  kinetic  energy  per 
unit  mass  as  a  function  of  depth.  At  the  thermocline  the  ratio  is  ap¬ 
proximately  3  to  1  while  at  1500  db  it  is  8  to  1. 

The  APE  per  unit  area  (or  TGPE)  between  300  and  1100  db  has  a  mean 

-4  -2 

value  over  the  eleven  cruises  of  .5  x  10  J  m  ,  with  a  standard  de~ 

-4  -2 

viation  of  .2  *  10  J  m  (Table  2.5).  Bray  and  Fofonoff  (1980)  calcu¬ 
lated  TGPE  between  300  and  2500  db  for  five  two-week  time  windows  during 

-4  -2 

MODE;  the  mean  value  was  .8  x  io  J  m  with  standard  deviation  of 

-4  -2 

.1  x  10  J  m  . 

The  change  in  potential  energy  of  the  system  corresponding  to  changes 
in  the  reference  steric  field  from  one  cruise  to  the  next  was  also  calcu¬ 
lated  (Table  2.5).  That  change  in  potential  energy  shows  no  trend  over 
the  three  year  period. 

In  Chapter  3  it  was  demonstrated  that  the  observed  seasonal  signals 
in  salinity  anomaly  and  vertical  displacement  are  consistent  with  a  model 
which  hypothesizes  forcing  by  large  scale  periodic  windstress  patterns  on 
a  3-plane  ocean  with  a  surface  mixed  layer,  depth  variable  stratification 
and  a  meridional  eastern  barrier.  The  historical  meteorological  data  for 
the  North  Atlantic  was  used  to  estimate  the  forcing  parameters,  and  the 
stratification  and  mixed  layer  parameters  were  determined  from  the  hydro- 
graphic  data.  The  model  ocean  response  to  the  atmospheric  forcing  has 
two  parts:  a  forced  wave  whose  horizontal  scales  are  those  of  the  forcing, 
and  free  waves  whose  meridional  scales  correspond  to  that  of  the  forcing 
and  whose  zonal  scales  are  determined  primarily  by  the  stratification. 

Only  the  barotropic  and  first  baroclinic  modes  are  not  coastally  trapped. 
The  long  baroclinic  mode  (wavelength  390  km)  is  found  to  be  the  only  mode 


if 
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which  can  be  important  in  the  interior  since  the  long  barotropic  mode  is 
very  much  longer  than  the  basin,  and  therefore  has  negligible  velocities 
and  vertical  displacements.  The  model  predicts  the  observed  phase  and 
horizontal  and  vertical  structure  of  both  salinity  anomaly  and  vertical 
displacements  within  the  errors  of  the  observations.  The  amplitude  of 
the  windstress  used  in  the  model  was  chosen  such  that  the  model  predic¬ 
tion  of  the  salinity  anomaly  signal  agreed  reasonably  well  with  the  ob¬ 
servations.  The  predicted  displacement  amplitude  agrees  with  the 
observations,  which  is  an  independent  test  of  the  consistency  of  the 
model.  The  model  further  predicts  changes  in  horizontal  and  vertical 
velocities,  APE  and  density.  The  predicted  APE  per  unit  mass  agrees  in 

amplitude  and  vertical  structure  with  the  annual  component  of  the  ob- 

2  2 

served  APE  per  unit  mass.  That  component  is  calculated  as  .  5N  a^  , 
with  a^  the  amplitude  of  the  seasonal  changes  in  vertical  displacement  it. 
The  changes  predicted  by  the  model  for  the  large  scale  density  field 
(corresponding  to  the  forced  wave)  are  a  factor  of  three  smaller  than  the 
observed  change  in  the  reference  steric  field,  and  have  different  phase 
and  vertical  structure.  The  predicted  meridional  velocities  are  a  maxi¬ 
mum  of  .018  m  sec  ^  at  the  base  of  the  mixed  layer  decreasing  to  .010  m 

sec  1  at  1000  db;  the  vertical  structure  is  that  of  the  first  baroclinic 

-5  -1 

mode.  The  zonal  velocities  predicted  by  the  model  are  of  order  10  msec 
With  regard  to  the  corollary  objectives  of  this  work,  it  was  found 
that  the  adiabatic  leveling  technique  described  by  Bray  and  Fofonoff 
(1980)  permits  the  precise  identification  of  surfaces  appropriate  for  ex¬ 
amining  tracer  distributions.  In  this  type  of  analysis,  it  is  important 
that  a  minimum  of  vertical  averaging  be  done  in  determining  potential 
density  surfaces  in  order  that  resolution  is  not  impaired,  but  sufficient 
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that  statistical  reliability  is  achieved.  The  analysis  of  observations 
presented  above  could  not  have  been  accomplished  using  "surfaces"  of  con¬ 
stant  potential  density  calculated  as  the  average  over  r  )mo  range  in 
potential  density:  too  much  information  is  lost  in  ordei.  to  reduce  the 
finestructure  errors.  The  adiabatic  leveling  technique  uses  a  least 
squares  method  to  weigh  the  information  available  from  all  data  points 
within  a  vertical  interval  to  determine  the  reference  steric  surface. 
Information  about  the  structure  of  the  density  field  is  thereby  retained, 
allowing  better  resolution  of  the  reference  surface  with  smaller  fine- 
structure  errors.  In  addition,  the  technique  provides  an  objective  esti¬ 
mate  of  those  errors.  An  important  further  advantage  of  the  adiabatic 
leveling  technique  from  the  perspective  of  data  analysis  is  that  vertical 
displacements  and  APE  are  also  calculated  accurately  using  the  reference 
steric  field;  the  technique  thereby  provides  a  consistent  overall  frame¬ 
work  for  the  simultaneous  analysis  of  observations  in  terms  of  tracers 
and  energetics. 

Given  the  distribution  of  some  tracer  along  the  reference  steric  sur¬ 
face,  it  becomes  necessary  to  develop  a  method  for  quantifying  that 
distribution — the  second  corollary  objective.  The  method  employed  in  this 
work  exploits  minimization  by  least  squares:  the  general  approach  is  to 
choose  a  possible  structure  or  sum  of  structures  (polynomials  or  sinu¬ 
soids)  and  calculate  the  amplitude  of  each  structure  using  a  least  squares 
fit  of  the  data.  This  approach  allows  the  determination  of  the  statistical 
confidence  that  the  chosen  structure  in  fact  exists  in  the  observations, 
and  what  the  error  in  the  determined  amplitude  is.  The  drawback  is  that 
the  bias  of  the  type  of  structure  expected  must  be  supplied  by  the 
investigator. 
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The  third  corollary  objective  was  the  accurate  estimation  of  APE, 
particularly  the  determination  of  the  errors  involv'd  in  that  estimate. 
It  was  found  that  the  Boussinesq  approximation  to  APE  per  unit  mass  is 
accurate  to  about  15%  in  most  applications;  the  error  introduced  by  neg¬ 
lecting  the  effects  of  compressibility  is  not  generally  significant. 
However,  errors  in  the  numerical  integration  over  pressure  in  the  cal¬ 
culation  of  APE  can  be  large  if  the  vertical  structure  is  not  resolved, 
so  that  care  must  be  taken  in  choosing  the  vertical  separation  of  points 
used  in  the  integration.  In  cases  where  the  vertical  displacements  of 
density  surfaces  are  large,  the  Boussinesq  approximation  may  become  in¬ 
valid;  there  the  gravitational  potential  energy  (GPE)  calculation,  which 
is  exact,  must  be  used  in  preference  to  the  Boussinesq  which  is  one  term 
in  a  Taylor  series  expansion  about  displacement.  Care  must  also  be  ex¬ 
ercised  in  the  definition  of  the  reference  steric  field:  if  the  ref¬ 
erence  field  used  does  not  conserve  mass  in  the  transformation  to  the 
observed  field,  large  systematic  errors  in  the  calculation  of  GPE  may 
result.  Thus,  in  calculating  GPE  for  an  isolated  feature  i.uch  as  a  Gulf 
Stream  ring,  the  density  field  away  from  the  ring  cannot  be  used  as  the 
reference  field  (Reid,  Elliott,  and  Olson,  1980) . 


4.2  Implications 

The  steady  circulation  of  the  ocean  is  presumed  to  be  driven  pri¬ 
marily  by  the  steady  component  of  the  windstross.  The  seasonal  fluctu¬ 
ations  in  the  windstress  are  of  the  same  order  of  magnitude  as  the  mean 
annual  values,  which  suggests  that  there  may  be  a  significant  seasonal 
component  to  the  general  circulation.  The  observational  and  theoretical 
results  described  in  the  previous  section  suggest  that  there  are  seasonal 
changes  of  meridional  velocity  in  the  eastern  North  Atlantic  of  order 
+.01  m  sec  *  with  zonal  wavelength  of  about  400  km  and  meridional  wave¬ 
length  of  perhaps  6000  km.  The  short  zonal  scale  suggests  that  horizon¬ 
tal  averaging  could  replace  temporal  averaging  in  distinguishing  this 
annual  signal  from  the  long-term  mean  flow.  RecalJ  from  Chapter  1  that 
the  mean  interior  flow  in  the  outer  Bay  of  Biscay  region  has  been  esti¬ 
mated  as  less  than  .01  m  sec  1  by  Swallow,  et  al.  (1977). 

The  relationship  of  horizontal  scale  to  amplitude  of  measurable  quan¬ 
tities  such  as  displacements  and  velocities  is  central  to  this  observa¬ 
tional  work  and  its  interpretation.  A  finite  total  energy  is  transferred 
from  the  atmosphere  to  the  ocean  during  the  seasonal  cycle.  That  energy 
is  transferred  over  large  horizontal  scales,  so  that  the  associated  dis¬ 
placements  are  relatively  small.  If,  by  some  unspecified  mechanism,  all 
or  part  of  that  energy  is  converted  to  smaller  horizontal  scales,  the 
relative  magnitude  of  the  displacements  will  increase.  The  central  hypoth¬ 
esis  of  this  work  is  that  just  such  a  mechanism  is  found  in  the  reflection 
of  a  large  scale,  atmospherically  forced  wave  at  a  meridional  barrier. 

The  shortness  of  horizontal  scale  wnich  permits  the  observation  of  these 
seasonal  signals  limits  their  direct  effect  upon  large  scale  oceanic 
phenomena  such  as  the  general  circulation  and  the  zonally  averaged 
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meridional  heat  flux.  However,  the  smaller  scale  annual  signals  can  in¬ 
fluence  the  interpretation  of  measurements  of  large  scale  quantities  un¬ 
less  adequate  spatial  or  temporal  averages  are  used.  For  that  reason 
alone  it  is  important  to  know  the  expected  scales  and  corresponding  amp¬ 
litudes  of  the  annual  signal  in  the  ocean. 

The  study  of  large  scale  circulation  through  examination  of  tracer 
distributions  has  a  long  history  in  oceanography:  initially,  because 
direct  current  meter  measurements  were  not  available,  and  more  recently 
because  the  large  scale  tracer  distributions  represent  a  time  integral  of 
the  velocity  field,  with  the  small  scale  variability  in  velocity  smoothed. 
With  the  observations  now  available  in  which  tracer  distributions  can  be 
studied  as  a  function  of  time,  the  low  frequency  variability  of  the  gen¬ 
eral  circulation  can  be  inferred  by  techniques  like  those  used  in  this 
work.  The  method  used  here  is  appropriate  to  regions  in  which  there  are 
steady  horizontal  gradients  of  tracers;  there  are  a  wide  range  of  condi¬ 
tions  under  which  the  method  may  be  usefully  applied.  Figure  4.1  is  an 
illustration  of  possible  solutions  of  the  one  dimensional  equation  for 
the  conservation  of  salinity  anomaly: 


S  +  v(x) S  =  0 
ti  y 


with  different  values  of  v  plotted  as  functions  of  and  S^.  (Values  of 
and  Sy  which  are  less  than  the  probable  errors  at  frequency  1  cpy  and 
scale  500  km  are  hatched.)  For  example,  in  order  to  detect  by  this  in¬ 
direct  a  method  a  velocity  signal  of  amplitude  1  cm  sec  and  frequency 


1  cpy,  a  steady  gradient  of  at  least  .2  ppm  km  is  necessary. 


Finally,  the  calculation  of  APE  presented  in  this  work  has  some  im¬ 
plications  for  the  study  of  ocean  energetics.  Previous  work  using  hydro- 
graphic  data  from  the  MODE  experiment  (Kim  1975;  Bray  and  Fofonoff,  1980) 


Figure  4.1  Particle  velocities  as  a  function  of  lateral  salinity 

anomaly  gradient  and  change  of  salinity  anomaly  with  time 
from  equation  (4.1).  Hatched  areas  represent  values 
less  than  the  probable  errors  for  frequency  1  cpy  and 
horizontal  scale  500  km. 
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indicated  a  simple  vertical  structure  for  APE  per  unit  mass,  resembling 
the  average  buoyancy  frequency  profile,  which  has  a  maximum  at  the 
thermocline,  and  decreases  monotonically  below.  The  amplitude  of  the 
MODE  APE,  which  was  calculated  over  a  region  400  x  400  km  in  size,  no¬ 
where  in  the  water  column  exceeded  the  observed  low  frequency  kinetic 
energy  (KE)  by  more  than  a  factor  of  3.  In  the  Bay  of  Biscay  observa¬ 
tions,  APE  per  unit  mass  has  a  second  maximum  below  the  thermocline  of 
approximately  equal  amplitude;  the  amplitude  and  the  deep  vertical  struc¬ 
ture  both  change  significantly  with  time.  Furthermore,  the  ratio  of  APE 
to  KE  below  the  thermocline  in  the  Bay  of  Biscay  observations  (Fig.  2.6) 
is  larger  than  for  MODE,  with  a  value  at  1200  db  of  8:1.  The  differ¬ 
ence  may  be  due  in  part  to  the  larger  area  over  which  APE  is  calculated 
in  the  Bay  of  Biscay.  Nevertheless,  it  appears  that  there  is  an  excess 
of  APE  over  KE  in  this  part  of  the  eastern  North  Atlantic,  as  well  as 
significant  low  frequency  variations  in  the  vertical  structure  and  amp¬ 
litude  of  APE  below  the  thermocline. 
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4.3  Further  Work 

The  observations  presently  available  limit  further  progress  on  the 
problem  of  the  oceanic  response  to  seasonal  atmospheric  forcing:  re¬ 
finement  of  the  theoretical  work  presented  in  this  thesis  is  limited  in 
its  usefulness  by  the  sophistication  and  completeness  of  the  observa¬ 
tions  to  which  the  theory  is  applied.  An  array  of  hydrographic  stations 
with  more  dense  spatial  coverage  might  allow  comparison  with  more  compli¬ 
cated  models  (involving,  for  example,  the  coastline  geometry  in  more 
detail) ;  a  shorter  sampling  period  would  permit  the  study  of.  higher  fre¬ 
quency  phenomena;  a  longer  time  series  would  give  not  only  more  statis¬ 
tical  confidence  to  the  result,  but  would  provide  evidence  of  interannual 
variability  of  the  signal.  The  meteorological  observations  presently 
available  should  be  examined  for  details  of  the  structure  of  the  low  fre¬ 
quency  wind  field,  with  particular  attention  to  the  direction  of  propaga¬ 
tion.  An  update  of  the  meteorological  observations  between  1972  and  the 
present  would  be  a  valuable  addendum  to  the  work  of  Bunker  and  Worthington 
(1976) ,  especially  for  synoptic  comparison  of  the  wind  field  and  the  large 
scale  oceanographic  experiments  which  have  been  carried  out  in  that 
interval. 

The  general  technique  for  data  analysis  presented  here  might  be 
applied  with  success  to  the  equatorial  oceans,  where  the  expected  signal 
is  strong;  as  an  example,  in  the  western  equatorial  Indian  Ocean  there 
are  strong  meridional  gradients  of  salinity  anomaly  resulting  from  the 
intrusion  of  high  salinity  water  from  the  Arabian  Sea.  Hydrographic  and 
velocity  measurements  of  considerable  extent  and  duration  are  available 
for  that  area  (Luyten  and  Swallow,  1976) ,  and  the  theory  of  low  frequency 
equatorial  waves  is  well  developed  (Philander,  1978) .  Eriksen  (1980) 
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discusses  observational  evidence  of  equatorially  trapped  waves  in  the 
western  part  of  the  Indian  Ocean,  using  moored  current  meter  data. 

The  adiabatic  leveling  technique  would  be  improved  by  extension  to 
depths  shallower  than  300  m,  where  it  is  presently  invalid.  Such  an 
extension  requires  that  a  method  be  developed  for  handling  density 
surfaces  which  are  not  continuous  over  the  region  considered.  A 
similar  problem  may  be  encountered  at  depth  in  regions  of  abrupt,  iso¬ 
lated  topographic  features. 
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Appendix  A:  Calculation  of  Systematic  and  Random  Error  in  APE 

1.  Systematic  differences  between  GPE  and  APE^ 

The  inclusion  of  compressibility  effects  in  GPE  was  shown  in 

Chapter  1  to  introduce  systematic  differences  between  GPE  and  APER. 

Those  differences  were  interpreted  as  contributions  from  horizontal  and 

vertical  gradients  of  compressibility,  and  it  was  noted  that  -a*?  can 

P 

be  replaced  by  kit  in  the  horize  -.tal  term.  It  should  also  be  noted  that 
the  computational  algorithm  calculates  ap(Pj^  which  differs  from 
a*(Pf).  Differentiating  (1.14)  with*^  =  a: 


da. 
_ i 

dp 


da„ 

(P.)  =  -rf 
l  dp 


«*£>  11  +  I’ 


+  K(pJ  -  K(pf)  [1  + 


dTT. 

dpJ 


or 


a*  (p. )  =  a*  (p,)  (1  + 

p  i  p  f 


d7T. 
dp  ' 


Thus  a*ir  may  be  rewritten  as 
P 


a*(p  )ir  =  a*(p.)ir  + 

p  f  pi 


_d_ 

dp 


*r 

J2_ 
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IT2 
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da* 

_R 

dp 


(A  1) 


This  expansion  may  also  be  used  to  evaluate  the  contribution  from  hori¬ 
zontal  gradients  of  compressibility;  however,  errors  in  the  determination 
of  the  second  derivative  of  a  make  it  less  satisfactory  than  the  simpler 
ktt.  The  corresponding  errors  in  TGPE  are  given  by 
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horizontal  gradients 


P  P 
r  'r 

■I  if. 

P  P 


dp  dp' 


P  P  - 

r*  Xr  p  „2 


vertical  gradients: 


-J  ?/  -i-*"*' 


P  P 


The  contributions  of  these  terms  to  TGPE  between  300  and  1100  db  are 
given  in  Table  2.5;  the  contributions  to  GPE  for  Phygas  42  are  plotted 
in  Pig.  2.6. 

2 .  Random  errors 

A  general  expression  for  the  variance  of  HPE  (1.11)  may  be  derived 

by  assuming  that  all  errors  are  uncorrelated  from  station  to  station  and 

from  level  to  level  in  the  vertical,  and  that  the  variance  of  a*  is 

P 

small  compared  to  that  of  the  displacement  ir^.  Then,  in  terms  of 
sums  over  the  integration  steps  Ap,  the  variance  of  HPE  is: 

j  'V  M 


V(ITO  -  vt  X  [5  Y.  2g  %  "k1  .  4  V 


j=j (p)  k=l 


j(P  ) 
r 


M  2 
_  a** 


-I  £1  Jrv(<”  4pr 


(A  2) 


j=j (p)  k=l 


Here  M  is  the  number  of  stations  and  the  subscript  j  refers  to  the 
level.  The  variance  of  it2  is  given  by 


V(TT2)  =  2[V(7Tk)]2  +  4  <TTk>2  Vf^) 


(A  3) 
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»  with  <ir,  >  the  ensemble  mean  of  tt,  ,  which  is  estimated  by  tt,  itself. 

k  k  k 

The  rmestructure  variance  of  IT,  is  taken  to  be  the  variance  of  the 
■  k 

regression  estimate  p,  which  is  described  in  Chapter  2.  The  finestruc- 
ture  errors  in  HPE  are  found  in  Table  2.5. 

Expressions  for  the  variance  of  TGPE  and  GPE  are  most  readily 
derived  using  the  expansion  of  (0^  -  d^)  from  equation  (1.14).  (This 
is  because  a  arid  af  are  obtained  by  inverting  high  order  polynomials.) 
If  any  covariance  between  the  terms  in  the  perturbation  expansion  of  TGPE 
are  ignored,  the  variance  of  TGPE  is  then  estimated  ass 


j(pr> 
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•  [mZt  v(<>  +  v£Z  °?)a#  Ap2j 

j= j (p)  k-l 
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(A  4) 

a 

with  V (tt  )  given  by  (A  3).  The  third  term  in  (A  4)  requires  some 
care  in  evaluation,  because  if  is  constrained  to  be  zero.  From  section 
1  of  the  Appendix: 

-Ct*7T  =  KTT. 

P 

a  is  very  nearly  constant;  if  it  is  taken  to  be  a  constant,  the  variance 
of  ktt  becomes 


V (KTT )  =  V(K  •  TT) 
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=  .K2  V(¥)  . 

In  practice,  7T  is  not  found  to  be  exactly  zero,  due  to  numerical  errors,  * 

so  that  some  small  contribution  is  to  be  expected  from  the  third  term. 

The  variance  of  it  is  estimated  as  (tt)2.  Values  of  the  finestructure 
errors  for  TGPE  are  found  in  Table  2.5  and  for  GPE  in  Fig.  2.6.  Because 
the  finestructure  errors  in  the  compressibility  terms  are  small,  the 
error  in  APE^  differs  very  little  from  that  in  GPE. 

1 

Measurement  errors  for  GPE  and  TGPE  were  calculated  by  assuming 
that  only  the  pressure  measurement  error  contributes  to  errors  in  tt 

K 

(i.e.,  that  temperature  and  salinity  errors  do  not  contribute)  and  that 

measurement  errors  in  a*,  and  F,  may  be  ignored.  A  value  for  the 

p  x 

pressure  error  of  ±  5  db  was  used  as  the  standard  deviation  of  tt  in 

A 

(A  5) .  Contributions  to  the  total  measurement  error  from  the  error  in 
the  compressibility  terms  were  found  to  be  negligible;  hence  the  measure¬ 
ment  error  in  APEB  is  effectively  the  same  as  that  calculated  for  GPE 
using  (A  4).  The  measurement  error  is  plotted  for  GPE  in  Fig.  2.6  and 
tabulated  for  TGPE  in  Table  2.5. 
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Appendix  B 

Cubic  Spline  Coefficients  for 

Iselin  - 

Worthington  -  Metcalf 

0-S  Curve 

(L.  Armi  and 

C.  Maillard,  personal 

communication) 
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Appendix  B  (cont. ) 


Salinity  is  then  calculated  as 


sn  =  c  +  c  At  +  c  At2  +  e  Atj 


with  AT  =  T  ,  -  T 

obs  0 
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